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ABSTRACT 


A  case  of  explosive  extratropical  cyclogenesis  in  the 
North  Pacific  is  analyzed  employing  quasi-Lagrangian  diag- 
nostic techniques  in  isobaric  coordinates.  The  First  GARP 
Global  Experiment  (FGGE)  Level  Ill-b  data  set  from  the  Euro- 
pean Center  Fcr  Medium- range  Heather  Forecasting  (ECMWF)  is 
used  in  a  synoptic  investigation  and  mass  and  vorticity 
budget   diagnostic  evaluation   of    storm   development. 

Mid-trcpospheric  positive  vorticity  advection  (PVA)  acts 
in  concert  with  low-level  thermal  advection  to  enhance 
surface  layer  organization  during  the  initial  periods. 
Explosive  development  occurs  under  upper -level  zonal  flow 
when  the  surface  system,  developing  along  the  low-level 
baroclinic  zone,  moves  under  the  front  left  quadrant  of  a 
strong      jet   streak.  Stability      in      the    lower      troposphere 

decreases  steadily  as  an  intense  two-layer  mass  circulation 
is    established.  The  divergence   term   provides      the   largest 

contribution  to  the  vorticity  tendency  during  explosive 
development  due  to  substantial  surface  layer  convergence  and 
upper-level   divergence. 
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I.        INTRODUCTION 

Extratropical  maritime  cylones  restrict  commercial, 
industrial   and   military      use   of    the   sea.  The    more    intense 

storms  which  form  over  very  small  time  periods  are  espe- 
cially dangercus  to  navigation  due  to  the  effects  of  the 
associated  wind,  waves  and  visibility  restrictions.  Accu- 
rate and  reliable  prediction  of  these  atmospheric  processes 
has  teen  a  goal  for  many  decades.  During  the  past  few  dec- 
ades, significant  advances  in  forecasting  these  synoptic- 
scale  atmospheric  circulation  features  have  been  realized 
through  the  development  of  numerical  weather  prediction 
(NWP)    models. 

While  current  numerical  models  provide  satisfactory 
accuracy  in  forecasting  most  synoptic-scale  processes,  they 
perfcrm  pocrly  in  predicting  the  predominantly  cold-season 
explosive  genesis  of  maritime  cyclones.  Sanders  and  Gyakum 
(1980)  define  explosive  cyclogenesis  as  an  extratrcpicai 
cyclone  which  undergoes  central  pressure  falls  of  at  least  1 
mb/h  for  24  h.  Accurate  forecasting  of  these  meteorological 
events  is  liiited  by  the  scarcity  of  observations  over 
oceanic  areas,  which  inhibits  the  model  initialization  pro- 
cess, and  by  the  smaller  spatial  and  temporal  scales  of  the 
explosively   developing  cyclone. 
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Accurate  forecasting  of  these  maritime  events  is  of 
great  significance  to  commercial,  industrial  and  military 
organizations  which  routinely  use  the  oceans  in  the  trans- 
action of  business  and  in  the  fulfillment  of  commitments. 
Several  recently  publicized  examples,  specifically,  the  less 
of  the  Russian  trawler  Metros  toy  (Mariners  Weather  Leg, 
1979)  ,  and  the  sinking  of  the  oil  rig  Ocean  Ranger  (Lemoyne, 
1982)  ,  illustrate  only  a  minute  fraction  of  the  total  less 
of  life  and  operational  hazards  associated  with  these 
intense    meteorological  events. 

Explosive  deepening  is  a  characteristic  of  the  vast 
majority  of  the  most  intense  cyclones  (Sanders  and  Gyakum, 
1980).  In  connection  with  associated  severe  weather  and  the 
failure  of  numerical  forecast  models  to  accurately  predict 
these  events,  it  is  of  paramount  importance  that  the 
processes  which  contribute  to  the  formation  of  these  vigor- 
ous storms  be  examined  and  fully  understood.  Several  recent 
studies  have  focussed  attention  on  detailed  analyses  of 
these  processes.  Roman  (1981)  completed  a  mass  budget  study 
and  Conant  (1982)  completed  a  circulation  and  angular  momen- 
tum budget  study  of  a  case  of  east-coast  cycloganesis. 
In tercomparisens   cf      these  results    with      conclusions    derived 
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in  this  thesis  provide  valuable  insight  into  the  formative 
mechanisms  of  maritime  explosive  cyclogenesis  and  lead  to  a 
general  quantitative  understanding  of  the  forecast  problem 
associated   with   these  intense  systems. 

This  thesis  is  part  cf  a  broadly-based  research  effort 
attempting  to  understand  better  the  formative  processes 
inherent  in  the  genesis  of  maritime  cyclones.  Specifically, 
this  study  concentrates  on  a  clearly  defined  structured 
approach  to  the  analysis  of  a  case  of  explosive  cyclogenesis 
over  the  North  Pacific  Ocean  during  the  period  12-17  January 
1979.  The  mass  and  circulation  (vorticity)  budgets  are 
examined  using  quasi-Lagr angian  diagnostics  (QLD)  as  devel- 
oped by  Jchnscn  and  Dcwney  (1975  a  and  b)  in  order  to  deter- 
mine the  role  cf  dynamic  forcing  in  the  budget  volume.  The 
local  contributicn  of  vorticity  sources  and  sinks  and  the 
vertical  and  horizontal  mass  and  vorticity  transports  are 
specifically   explored. 

Identification  of  these  processes  is  uniquely  interest- 
ing in  that  the  cyclcne  under  study  consists  of  the  classic 
frontal  wave  development  from  strong  low  level  baroclinic- 
ity.  The  distinguishing  feature  is  the  absence  of  a  strong 
migratory    short    wave    aloft,      which   categorizes   this   storm    as 
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a  Pettsrssen  type-A  system,  as  opposed  to  the  type-B  fea- 
ture which  contains  a  500  ib  short  wave  responsible  for  the 
rapid  surface  deepening  (Eetterssen  and  Smebye,  1971).  in 
the  type-A  case,  the  flow  aloft  is  predominantly  zonal 
through  the  early  explosive  deepening  phase  with  positive 
vorticity  advection  of  minor  significance.  Low  level  ther- 
mal advection  plays  the  important  role  in  -he  deepening 
process.  This  thesis  is  believed  to  be  one  of  the  first 
detailed    analyses     of  this     class  of      maritime   storms.  In 

addition,  the  vast  majority  of  explosive  deepeners  form  in 
the  North  Pacific  Ocean  (Sanders  and  Gyakum,  1980).  Inas- 
much as  the  most  severe  weather  is  directly  associated  with 
these  severe  storms,  this  study  assumes  particular  signifi- 
cance in  attempting  to  understand  this  meteorological 
phenomenon    in    crder    to    mirimize    its   disastrous    effects. 

The  specific  approach  employed  in  this  thesis  is  the  use 
of  QLD  techniques  in  an  isobaric  spherical  coordinate  system 
which   is      translated    with   the   stcrm.  Horizontal   advection 

associated  with  wave  translation  is  isolated  from  the  rela- 
tive advection  and  the  divergence  of  the  transport  associ- 
ated with  cyclone  development.  The  fluctuations  of 
meteorological   variables      are   then    analyzed    through      the   use 
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of    a   cylindrical     budget    volume    centered   on      and   translating 
with   the      cyclcne  surface   center.  Vertical   distributions, 

lateral  exchanges  and  sources  and  sinks  of  cyclone  proper- 
ties resulting  from  purely  developmental  processes  are  then 
analyzad.  Application  of  QLD  techniques  in  studying  cyclone 
development  (Wash,  1978)  has  proven  valuable  in  diagnoses  of 
poorly  forecast  storms,  and  Roman  (1981)  and  Conant  (1982) 
have  demonstrated  QLD  applicability  in  detailed  analyses  of 
east-coast   cyclogenesis. 

The   main    objectives    of  this    thesis  are: 

•      A  synoptic    investigation   of      the   classical    frontal    w< 
cyclogenesis  in    the    North   Pacific  Ocean; 
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•   A  diagnostic  investigation  of   storm  dynamics  through 
the  use  of  QLD  budget  analyses  of  mass  and  vorticity. 

The  detailed  analyses  of  these  processes  is  further 
aided  by  an  exceptionally  strong  supporting  data  base.  The 
case  under  study  occurred  during  the  first  special  observa- 
tion period  (SOP-1)  of  the  First  GARP  Global  Experiment 
(FGGE)  .  This  FGG2  data  set  was  assimilatad  by  the  European 
Center  for  Medium-range  Weather  Forecasting  (ECMWF)  and 
obtained  through  the  auspices  of  the  Goddard  Space  Flight 
Center.   This  data  set  is  special,  in  that  observations  from 
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aircraft,  ships,  buoys  and  satellites  gathered  during  an 
intense  collection  effort  were  manually  checked  and  then 
systematically  integrated  into  a  dense  network  of 
observations.  Consequently,  the  supporting  data  base  ana- 
lyzed in  this  thesis  effort  is  exceptionally  broad-based, 
complete,  and        integrated  into        a         well-diversified 

observational   network. 

Quantitative  analysis  of  FGGE-type  data  bases  gives  spe- 
cific insights  into  development  processes  inherent  in  mari- 
time explosive  cyclogenesis.  In  the  following  chapters, 
examination  of  the  FGGE  data  (Chapter  2)  and  an  abbreviated 
description  cf  QLD  procedures  (Chapter  3)  are  presented. 
Chapter  4  discusses  the  cyclone  under  study  through  a  synop- 
tic overview  which  highlights  significant  flow  patterns  and 
associated  upper  and  lower-level  features.  Chapters  5  and  6 
present  the  results  of  the  budget  calculations  and  Chapter  7 
contains  specific  conclusicns  from  these  budget  calculations 
and    discusses   recommendations   for   further    research. 
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II.       ENVIRONMENTAL    DATA    BASE    DESCRIPTION 

A.        FGGE    DATA 

Curing  the  first  special  observation  period  (SOP-1)  of 
the  Global  Weather  Experiment,  the  largest  concentration  of 
meteorological  resources  ever  assembled  was  deployed  to 
observe  the  atmosphere  and  oceans.  The  quality  of  the  data 
generated  from  this  effort  marked  a  milestone  in  the 
development  of  numerical  weather  prediction  models  and 
synoptic    meteorology      in    general.  By  any      standard,       this 

experiment  shewed  an  unparalleled  example  of  scientific 
cooperation   and   excellent   resource    management. 

The  FGGE  Level  Ill-b  data  set  used  in  this  thesis  was 
constructed  frcm  observations  integrated  through  a  data 
assimilation   system      at    ECMWF.  Bengsston    et      al.,       (1982) 

describe  the  system  as  a  three-dimensional,  multivariate, 
optimum  interpolation  scheme  using  a  nonlinear  normal-mode 
initialization  and  an  automatic  system  for  data  checking. 
The  direct  field-level  observations  (Level  I-b)  were  trans- 
formed into  the  basic  meteorological  variables  (Level  Il-fc) , 
which  constituted     the    basic     meteorological    product      of    the 
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experiment.  A  15-ievel  model  with  horizontal  resolution  of 
1.875  degrees  of  latitude  and  longitude  was  used  to  trans- 
form the  level  Il-b  data  to  Level  Ill-b.  At  the  Ill-b 
level,  the  observations  were  manually  checked  and  objec- 
tively analyzed  into  a  set  of  dynamically  derived  meteorolo- 
gical fields.  The  flowchart  of  this  data  assimilation 
process  is  depicted  in  Fig.  1.  The  FGGS  Level  Ill-b  data 
production  effort  was  completed  in  June  1981  and  the  data 
fields  were  catalogued  for  6-h  observations  for  all  standard 
levels  up  to  and  including  10  mb.  A  detailed  description  of 
this  assimilation  scheme  can  be  found  in  3engsston  at  al., 
(1982)  . 

The  FGGE  data  base  is  currently  undergoing  validation  by 
a  variety  cf  groups.  Roman  (1981)  and  Conant  (1982)  demon- 
strated its  validity  on  a  case  of  east-coast  cyclogenesis. 
In  an  independent  assessment  of  the  FGGS  data  set,  Halem  et 
al.,  (1982)  present  a  detailed  analysis  in  connection  with 
the  coverage  and  accuracy  of  the  remote  sensing  portion  of 
the  network  covering  the  oceanic  regions.  Analyses  based  en 
these   inputs    provided  the    following  conclusions: 
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•  The  FGGE  ctserving  system  has  a  positive  effect  on 
forecast   skill    using   Si    criteria; 

•  Significant  corrections  to  the  300  mb  first  guess 
fields   occur  over   oceanic   regions   with   the   FGGE   data; 

•  Nearly  all  of  the  forecasts  that  retain  useful  skill  by 
day    five   shew  a    positive    impact   from    the   FGGE  system; 

•  In  certain  synoptic  situations,  the  FGGE  observina  sys- 
tem is  capable  of  measuring  the  atmosphere  alobaily 
with  sufficient  accuracy  for  current  NWP  models'to  make 
useful   extended- range  forecasts   beyond  one    week. 

These  and  other  conclusions  as  well  as  a  detailed  descrip- 
tion of  the  analysis  scheme  can  be  found  in  Halem  et  al., 
(1981)  and  Baker  et  al. ,  (1981)  .  It  is  evident  from  these 
and  ether  studies  that  the  FGGE  Level  Ill-b  data  set  pro- 
vides a  strong  foundation  for  accomplishing  the  objectives 
of    this   thesis. 

A  summary  of  the  FGGE  data  availability  by  field-type 
and  level  is  depicted  in  Table  I.  The  advantages  of  this 
type  of  augmented  data  support  over  the  standard  data  base 
include: 

•  A  considerably  richer  data  base,  particularly  over 
data-sparse   oceanic    areas; 

•  The  high  resolution  grid  is  convenient  for  deoicting 
the  smaller  synoptic-scale  features  associated  with 
cyclcgenesis; 

•  Easic  and  derived  fields  (vertical  velocity,  humidity, 
etc.)       are    included    in  the    data    set; 
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•  The      full   area      of    cyclogenesis      and    storm      translation 
enjoys   si  Hilar    data    coverage. 

These  significant     advantages    demonstrate      the   value      of    the 

FGGE   Level   Ill-b   data   set      in   quantitative    analyses   of    storm 

features. 

B.        SUMMARY    OF    ANALYSES    AND    PROGNOSES    UTILIZED 

Three    sgzs   of  analyses  are   utilized  in    this   thesis: 

•  National    Meteorological   Center    (NMC)     analyses   and    prog- 
noses; 

•  Fleet    Numerical    Oceanography  Center    (FNOC)    analyses   and 
prognoses; 

•  ECMWF    Level    Ill-b   analyses. 

These  analyses  are  used  in  conjunction  with  the  synoptic 
discussion  (Chapter  4)  and  for  the  verification  of  the 
forecast  problem  alluded  to  earlier.  A  complete  listing  of 
all  analyses  and  prognoses  utilized  in  this  study  is  con- 
tained in  Table  II.  A  complete  description  of  the  FNOC 
products  can  be  found  in  the  U.S.  Naval  Weather  Service 
Numerical    Environmental    Products    Manual    (NAVAIR   50-1G-522) . 
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III.       C.UASI-LAGBANGIAN    DIAGNOSTIC     (OLD) 

Tggnrrsjns — 


To  derive  generalized  budget  analysis  techniques  for 
processes  related  to  cyclone  development,  changes  in  -he 
budget  volume  properties  due  to  migration  of  the  system  must 
be  separated  to  isolate  true  developmental  trends.  In  the 
past,  syncptic-scale  diagnostic  approaches  (e.g.  Rung  and 
Baker,  1975)  evaluated  storm  processes  relative  to  a  fixed 
coordinate  system.  This  type  of  approach  is  effective  for 
selected  time  periods  when  storm  translation  is  negligible, 
but  suffers  markedly  when  applied  to  rapidly  translating 
systems   of   appreciable   intensity.  In  particular,      changes 

may  cccur  due  tc  system  ccvement,  system  development,  or 
both.  The  only  accurate  approach  to  assess  true  developmen- 
tal fluctuations  associated  with  cyclone  development  is  tc 
construct  a  moving  budget  volume  centered  on  some  distin- 
guishing storm  feature.  In  this  way,  horizontal  advection 
associated  with  wave  translation  may  be  isolated.  In  addi- 
tion, any  vertical  redistribution  of  a  budget  property  may 
be  analyzed  in  terms  of  vertical  transport  processes  within 
the    cyclone    budget   regime.        QLD    techniques    employed    in    this 
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thesis  provide  the  means  fcr  this  type  of  approach.  Roman 
(1981)  and  Conant  (1982)  have  demonstrated  the  validity  of 
this   type    cf   approach   on    a  case    of   east-coast   cyclogenesis. 

To  study  the  fluctuations  in  any  storm  property,  a 
generalized  budget  equation  is  constructed  in  terms  of  the 
storm  budget  coordinate  system.  The  aquation  consists  of 
the  time  rate  of  change  of  the  desired  budget  quantity 
(df/dt)  and  the  lateral  flux  transport,  vertical  flux 
transport  and  sources/sinks  (see  Table  III) .  The  equation 
is  applicable  to  the  storm  budget  volume  in  spherical  coor- 
dinates. The  spherical  portion  of  the  atmosphere  incorpo- 
rated within  the  budget  volume  is  approximated  by  a  cylinder 
fcr  ease  in  computations  and  definition  of  budget  proper- 
ties. This  approximation  is  valid  for  the  relatively  thin 
slab  of  atmosphere  (1000-100  mb)  selected  for  budget 
analyses. 

The  formulation  of  the  storm  budget  coordinate  system 
and  details  cf  the  budget  equation  notation  are  discussed  by 
Wash  (1978).  The  storm  budget  volume  is  defined  in  terms  of 
10  mandatory  pressure  levels  in  the  vertical,  and  in  the 
horizontal  by  outwardly  directed  radii  originating  from  a 
central    point.         For    example,    the   budget    property   at   a    point 
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may  be  analyzed  by  describing  the  level/layer  (1000/925  jib)  , 
the  budget  radius  in  degrees  of  latitude  from  the  storm 
center  (6) ,  and  the  quadrant  desired  (grid  points  1-36) . 
The  mass  and  vorticity  budget  equations  employed  in  this 
analysis    are   depicted   in    Tables    IV   and     V    respectively. 

In  the  mass  budget  (Table  IV) ,  the  general  equation  rep- 
resenting the  change  in  mass  with  time  is  composed  of  con- 
tributions frcm  lateral  and  vertical  transports.  The 
lateral  transport  term  represents  the  net  horizontal  conver- 
gence of  mass  flux  within  the  budget  volume.  The  vertical 
transport  is  calculated  through  the  use  of  derived  vertical 
velocity  fields  and  is  cecessarily  a  reflection  of  net 
lateral  flux  into  the  volume.  Application  cf  boundary  con- 
ditions that  vertical  motion  is  zero  at  the  top  and  bottom 
surfaces  of  the  budget  vclume  simplifies  the  equation,  as 
shewn   in   the    middle    of  Table   IV. 

Several  useful  relationships  can  be  derived  frcm  this 
simplified  equation  and  the  boundary  conditions.  Contribu- 
tions to  the  mass  time  tendency  in  iso baric  coordinates 
arises  chiefly  from  pressure  fluctuations  in  the  surface 
layer,  as  the  other  layers  in  the  budget  are  fixed  pressure 
surfaces.  This      term        is      computed        easily     and        vith 
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considerable  accuracy.  The  net  lateral  mass  transports  in 
the  layers  are  calculated  from  the  wind  fields  and  necessar- 
ily balance  the  vertical  transports  through  continuity. 
Errors  in  the  wind  fields  cause  the  vertically  integrated 
lateral  mass  transport  to  differ  from  the  vertically 
integrated  mass  tendency,  and  cause  a  residual  in  the 
budget. 

Application  of  the  boundary  conditions  to  the  divergence 
of  the  vertical  transport  causes  the  vertical  sum  to  vanish, 
and  thereby  isolates  the  residual  correction  which  must  be 
applied  layer  by  layer  tc  force  the  vertical  sum  of  the 
divergence  tc  zero.  The  O'Brien  (1970)  linear  correction 
scheme  is  utilized  to  correct  for  these  residual  imbalances. 
The  upper  layers  are  weighted  more  due  to  the  greater  inac- 
curacies of  the  wind  measurements  at  these  levels.  These 
corrections  are  applied  tc  the  layer  mass  transports,  and 
are  combined  with  the  translation  correction  (discussed 
earlier)  to  give  the  final  corrected  form  of  the  mass 
transport    term.  The   vertical    velocities      representing   the 

vertical  mass  transport  are  then  recovered  through 
continuity. 
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In  the  vorticity  budget  formulation,  the  contributions 
to  vorticity  changes  occur  from  lateral  and  vertical  fluxes, 
sources      and      sinks       (Table      V,  Section      A) .  The      flux 

(transport)  terms  can  be  expanded  into  divergence  and  advec- 
tion  contributions  by  the  use  of  vector  identities  (Table  V, 
Section  B)  ,  which  facilitates  fur-her  insight  into  the  phys- 
ical processes  they  represent.  The  source/sink  terms  arise 
from  net  generation/dissipation  of  vorticity  within  the 
budget  volume,  and  in  iscbaric  coordinates  are  comprised  of 
divergence,  tilting  and  friction  components.  Notice  that 
three  divergence  terms  arise  from  this  expansion  of  the 
general  vorticity  budget  equation.  The  horizontal  transport 
component  is  identical  to  the  source  divergence  term,  and 
the  horizcntal  and  vertical  transport  components  are  equal 
and  opposite.  Thus,  the  net  contribution  of  divergence  ir. 
the  budget  is  manifest  ir.  the  source  term.  A  specific 
breakdown  of  the  partitioned  vorticity  budget  equation  is 
presented    in   Section    B   of    labia    V. 

To  isolate  the  contributions  of  both  the  rotational  and 
irrotational  wind  components  to  the  transport  of  a  desired 
budget      quantity,  mean      and      eddy      modes         (analogous      to 

divergence   and      advection    terms    ir.    the      vorticity   equation)  , 


26 


are  derived  frcm  the  transport  terms  (Table  VI)  .  The  pur- 
pose of  this  technique  is  tc  isolate  the  changes  in  the  mean 
budget  properties  from  these  occurring  on  the  eddy  scale. 
Thus,  a  comprehensive  analysis  of  the  lateral  and  vertical 
flux  terms  can  be  made  by  separating  them  into  mean  and  eddy 
mode  contributions  and  comparing  them  with  the  divergence 
and  advection  terms  for  analysis  of  the  physical  processes 
they   represent.  The  relation      of   the      mean/eddy   mode      and 

divergence/advection  partitions  to  the  total  transport  is 
also   illustrated   in    Section  B   of    Table   V. 

The    budget   analysis      begins    with   the   choice      of  the    par-* 
ticular      property  to      be    analyzed.  The    next      step    in      the 

budget  calculations  consists  of  interpolating  the  u  and  v 
compenents  of  the  ?GGE  data  set  onto  the  cylindrical  budget 
volume  grid.  Next,  the  interpolated  winds  are  converted  to 
tangential  and  radial  compenents  relative  to  the  outer  vol- 
ume boundary  and  corrected  fcr  storm  motion.  In  this  form, 
the  data  are  ready  to  be  used  in  the  budget  equations 
discussed  above.  Derived  fields  (moisture,  vertical  veloc- 
ity, kinetic  energy,  etc.),  not  directly  used  in  budget  com- 
putations can  be  used  tc  supplement  insight  into  storm 
dynamics.  Interesting      comparisons    between      FGGE      derived 
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vertical   velocity  and   QLD   calculated    vertical   velocity,       for 
example,    serve   to  validate  the  QLD   technique   approach. 

The   QLB   approach    provides   insight   into    the    following    key 
properties: 

•  The    magnitude      and   significance   of  the      horizontal    mass 
transport   structure; 

•  The    magnitude      and   significance  of  the      horizontal    vor- 
ticity   transport  structure; 

•  The    vertical  transport  structure     of    mass   and   vorticity 
derived    from  vertical   motion   estimates; 

•  The    magnitude      and   nature  of      source/sink   contributions 
to   budget   analyses. 

The   next    chapter     discusses  the    development   sequence      of    the 

storm     through      a      synoptic      discussion      utilizing      synoptic 

fields,    analyses   and    forecasts,    and   satellite   imagery. 
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IV.   SYNOPTIC  OVERVIEW 

A.   GENERAL 

Oceanic  cyclcgenesis  mechanisms  can  in  part  ba  portrayed 
by  graphical  representation  of  pertinent  meteorological  and 
oceanic  fields.  This  chapter  explores  the  synoptic  condi- 
tions associated  with  a  specific  case  of  North  Pacific 
explosive      cyclcgenesis      during       12-17     January      1979.  An 

attempt  is  mad€  to  gain  insight  into  formative  processes  and 
possible  factors  related  to  the  serious  shortcomings  in  the 
model  repesentaticn  of  the  dynamics  associated  with  these 
vigorous  events.  Use  is  made  of  FGGE  data  analyses.  Defense 
Military  Satellite  Program  (DMSP)  imagery  and  FNOC  sea 
surface   temperature     (SST)      analyses.  Comparisons  are   made 

between  FGGE  ,  FNOC  and  NMC  sea  level  pressure  analyses  and 
prognoses  to  highlight  differences  in  forecast  deepening 
rates  and  analyzed  and  forecast  storm  positions.  Explosive 
deepening  criteria  were  adopted  from  Sanders  and  Gyakum 
(1980),  who  addressed  the  NMC  7-layer  PE  model  performance 
associated  with  "explosive"  prediction.  Predicted  deepening 
rates   for    the      cyclones    they   studied   wer=      consistently    less 
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than  1/3  of  the  observed  rates,  ar.d  North  Pacific  explosive 
deepening  east  of  the  dateline  was  virtually  missed  by  the 
model.  An  assessment  of  the  FNOC  5-layer  PE  model  perform- 
ance (Gyakum,  1981)  demonstrated  similar  results. 
Specifically,  the  model  captured  only  1/3  cf  the  observed 
central  pressure  fails  for  all  time  periods.  In  addition, 
the  most  accurate  of  all  the  forecasts  were  shown  to  be 
valid  between  24-36  h.  In  comparing  the  FGGE  analyses  with 
the  FNOC  and  NMC  analyses,  it  should  be  noted  that  the  FNOC 
analyses  and  prognoses  are  derived  from  the  5-layer  PE 
coarse  mesh  (381  km  hcrizcntal  resolution)  model.  In  con- 
trast, the  NMC  analyses  are  strictly  hand-analyzed  Northern 
Hemisphere  strips,  while  the  prognoses  are  derived  from  the 
7-player    PE    fine    mesh     (190.5   km  horizontal    resolution)    model. 

B.        MODEL    FREEICTION    PERFORMANCE 

The  storm  track  (at  12-h  increments)  is  presented  in 
Fig.  2a.  The  cyclone  develops  southeast  of  Japan  in  the 
vicinity  30N,  135E  and  transits  toward  a  final  position  off 
the    west    coast   of  Alaska.  Notice   the  significant  transla- 

tion  increase   after    0000    GMT      13    January.         The    5-day    evolu- 
tion    (0000   GMT    12  January    1979      to    0000   GMT    17    January    1979) 
in   central      sea-level   pressure,       as      analyzed    from      the    FGGE 
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data,  NMC  and  FNOC  analyses,  is  illustrated  in  Fig.  2b. 
Assuming  the  representativeness  of  the  FGGE  data,  both 
series  of  FNOC  and  NMC  analyses  fail  markedly  in  represent- 
ing the  intense  deepening  phase  (0000  GMT  13  to  0000  GMT  15 
January  1979) .  There  is  close  agreement  between  FGGE  and 
NMC      analyses   at     other   times.  It      appears  that      specific 

inputs  to  the  FGGE  data  set  (mentioned  earlier)  not 
available  to  FNOC  or  NMC  for  real-time  analyses  provide  a 
more      accurate   portrayal      cf    the      deepening   period.  Close 

inspection  of  the  observations  included  in  the  FGGE  analyses 
shows  that  the  deepening  phase  is  well  documented  at  all 
levels.  At  1200  GMT  1U  January  1979,  for  example,  as  many 
as  seven  sea-level  observations  from  land  and  sea  stations 
ar 3  available  in  the  north  and  east  quadrants  of  the  stcrm 
and  form  the  surface  data  netwcrk.  In  the  upper  levels, 
satellite-derived  temperature  and  height  anomalies,  and  250 
mb  aircraft  winds  and  temperatures  combine  to  form  a  dense 
observational  network.  In  the  NMC  analysis,  ship  reports  in 
the  area  were  conspicuously  absent,  as  the  analysis  appeared 
to  be  based  upon  land  stations  only.  The  FNOC  analysis, 
however,  contained  four  of  the  seven  FGGE  surface 
observations,         tut    did      not      contain      the   two      observations 
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closest  to  the  stcrm  center,  and  this  caused  a  misrepresen- 
tation of  the  surface  central  pressure.  Thus,  the  more 
accurate  representation  of  the  deepening  phase  in  the  FGGE 
analysis  is  directly  attributed  to  a  superior  data  base  and 
not    just    tc   the   benefit    of  a    "better"    firs-    guess. 

A  comparison  between  the  36-h  sea-level  pressure 
forecasts  is  shown  in  Fig.  3a.  All  curves  show  the  general 
deepening  trends,  but  the  FNOC  and  NMC  forecasts 
consistently  lag  the  FGGE  analyzed  pressure,  and  are  highly 
variable  in  predicting  the  magnitude  of  the  central  pressure 
falls.  From  these  displays,  it  is  evident  that  both  models 
fail  significantly  in  the  timely  and  accurate  prediction  of 
storm  intensification,  especially  during  the  explosive 
stage.  This  example  illustrates  further  the  forecast  prob- 
lems discussed  by  Sanders  and  Gyakum  (1980)  alluded  to 
earlier. 

The  location  error  in  36-h  storm  positions  for  NMC  and 
FNOC  predictions,  relative  to  the  FGGE  analyzed  position,  is 
presented  in  Fig.  3b.  It  is  clearly  seen  that  significant 
errors  exist  in  both  forecasts  at  various  times  curing  the 
deepening    phase.  Such    significant   errors   in      the  analyzed 

and    forecast    features  of    this    system    reflect    an    inability   of 
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current  operational  models  to  accurately  forecast  these 
events. 

C.        SYNOPTIC    DISCUSSION 

The  synoptic  analysis  approach  in  this  section  consists 
of  an  examination  of  the  basic  meteorological  fields  in 
three  displays  of  two  fields  each.  Display  A  contains  250 
mb  winds  and  isotachs,  B  presents  500  mb  heights  and  abso- 
lute vorticity,  and  C  shows  1000-500  mb  thickness  and  sea- 
level  pressure  patterns.  This  approach  generally  allows 
only  subjective  insights  into  features  of  interest  during 
the   storm      life   cycle.  The    quantitative      analysis   of      the 

dynamical  processes  is  presented  in  the  following  chapter. 
The  intent  of  this  section  is  to  trace  the  development  of 
the  system  thicugh  a  synoptic  discussion  of  major  features 
contributing   tc   overall    storm   development. 

Examination  of  atmospheric  conditions  at  0000  SHT  12 
January  1979  in  the  area  southeast  of  Japan  shows  marked 
cyclcgenetic  potential.  The  surface  picture  (Fig.  He)  is 
dominated  ty  a  broad  trough  (vicinity  30N,  140E)  extending 
northeast-southwest  through  the  center  of  a  checkerboard 
pattern  of  fairly  strong  lew  and  high  centers.  The  signifi- 
cant   features    include  the    Siberian    High  near    55N,    100S    (just 
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off  the  western  edge  of  the  chart)  which  is  very  strong  and 
stationary.  The  parent  lew  in  the  vicinity  of  53N,  170E  is 
now  quasi-stationary  and  filling  (950  mb)  .  A  large  maritime 
high  (1030  mb)  dominates  the  north  central  Pacific  (vicinity 
30N,  178W) ,  and  the  northern  fringe  of  Typhoon  Alice  is  jus- 
visible  near  2QN,  140E).  The  combination  of  flow  from  these 
systems  serves  to  funnel  cold,  polar  continental  air  south- 
ward towards  northward  flowing  warm,  moist  tropical  air, 
which     enhances      the    frontogenes is      potential.  A      general 

destabilizaticn  cf  the  polar  air  occurs  north  of  the  area  of 
strong  low  level  baroclinicity  shown  by  the  concentration  of 
thickness    contours   just    east   of    Japan. 

In  the  mid-troposphere  (Fig.  4b)  ,  two  distinctive  fea- 
tures are  evident.  The  flow  over  the  western  Pacific  is 
dominated  by  an  extensive  region  cf  low  heights  associated 
with  a  quasi-stationary  leng  wave  trough,  which  is  a  portion 
of  a  four-wave  planetary  pattern.  This  significant  feature 
accounts  for  the  steady  progression  of  mobile  short  waves 
and  associated  surface  features  through  the  region  daring 
the  preceding  week,  in  addition  to  supporting  the  steady 
flow   of      cold   air  off     the   Siberian   continent.  The    second 

major   feature      of  interest  is      the    predominantly      zonal    flew 
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over     tha    entire     incipient  region      with   areas      of   weak      and 
small  absolute   vcrticity    upstream. 

In  the  upper  troposphere  (Fig.  4a)  ,  the  presence  of  a 
broad  jet  (  >  50  m/sec)  aloft  is  a  further  manifestation  of 
the  fcaroclinicity  of  the  lower  troposphere.  As  yet,  signif- 
icant jet  streak  (  >  75  m/sec)  propagation  into  the  region 
remains  well  upstream,  and  therefore  is  not  a  factor  in  the 
low-level    development. 

Examination  of  the  visual  (Fig.  5)  and  infrared  (Fig.  6) 
satellite  imagery  valid  at  0041  GMT  12  January  1979  shows  a 
couplet  of  similar  low  level  cloud  masses  (labelled  A  and  B) 
southeast        of      Japan.  Lack        of      significant        vertical 

davelopment  is  evident  on  the  infrared  imagery,  in  contrast 
to  the  area  cf  enhanced  convection  (labelled  C)  associated 
with  an  area  of  positive  vcrticity  advection  (?VA)  upstream 
over  Korea  (Fig.  4b) .  Cloud  mass  B  is  the  precursor  cf  a 
weak  cyclone  which  moves  to  the  east  and  rapidly  dissipates. 
Cloud  mass  A  is  the  precursor  of  the  cyclone  analyzed  in 
this  study.  The  evolving  low-amplitude  "s-shape"  configura- 
tion with  a  distinct  edge  on  the  western  flank  is  indicative 
of  increasing  thermal  and  low  level  absolute  vorticity  gra- 
dients  with   time    (Weldon,     1977).       Although   fairly    small,    the 
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cloud  mass  is  beginning  tc  take  on  -the  familiar  "leaf-shape" 
which  Weldon  (1977)  characterizes  as  typical  of  incipient 
oceanic   cyclogenesis. 

During  the  next  24  h  ending  at  0000  GMT  13  January  1979, 
surface  development  continues  (vicinity  35N,  150E)  intc  a 
closed  center  (995  mb) .  The  incipient  system  (Fig.  7c) 
appears  as  a  broad  area  of  low  pressure,  and  the  development 
of  a  weak  thermal  ridge  near  the  surface  center  indicates 
frontogenesis  continues  during  the  period.  Ships  in  the 
vicinity  reported  seas  on  the  order  of  15  ft,  with  periods 
ranging  between  3-5  s,  which  indicates  a  duration  limited 
wave  spectrum.  Of  interest  is  the  continued  strengthening 
of  the  Siberian  high  (now  1080  m  b)  which  steadily  pumps  cold 
air  southward,  in  association  with  the  backside  circulation 
of  the  stationary  cold  low  to  the  east.  Destabilization  of 
the  surface  layer  continues  as  cold  air  streams  southward 
over  the  warm  ocean  surface,  and  a  polar  low  has  become 
organized  ever  the  southeast  tip  of  Japan  due  west  of  the 
incipient  system  (indicated  on  detailed  NMC  analyses,  not 
shown   here) . 

In      the    aid-troposphere      (Fig.         7b)  ,         the    500      mb      low 
anchored    off    the    Kamchatka   Peninsula   has    deepened.         A    shcrt 
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wave  with  an  associated  broad  absolute  vorticity  area  has 
formed  over  Japan,  approximately  800  nm  upstream  of  the 
incipient  sysxem.  A  smaller,  less  distinct  vorticity  area 
has    formed   directly    over    the   surface  system. 

Increased  vorticity  in  the  mid-troposphere  can  also  be 
linked  to  the  approach  of  a  jet  streak  from  the  west.  At 
this  stage,  however,  the  position  of  the  jet  streak  indi- 
cates only  minimal  influence  on  the  low-level  circulation, 
as  the  jet  cere  is  directly  superposed  in  the  vertical,  and 
no  clear  pattern  of  cyclonic  shear  or  diffluence  can  be  sees 
at    250   mb     (Fig.    7a)  . 

The  ir.crsased  vertical  development  over  the  period  can 
be  seen  in  the  DMSP  satellite  imagery  (Figs.  8  ani  9)  valid 
at  2341  GMT  12  January  1979.  The  developing  system  is 
labeled  A  and  is  centered  approximately  600  nm  southeast  of 
Japan.  Ccld  plumes  extend  for  virtually  the  entire  breadth 
in  the  infrared  image  (Fig.  9)  ,  with  cloud  tops  becoming 
progressively  warmer  toward  the  tapered  end  to  the  south- 
west. Enhanced  convecticn  and  continued  development  of  the 
leaf  structure  is  shown  in  the  visual  image  (Fig.  8) .  The 
region  under  the  leaf  appears  to  contain  marked  baroclinic- 
ity      in   a      zone    in      advance     of    the      propagating    jet      streak 
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aloft,  which  is  similar  tc  the  conditions  observed  by  U el don 
(1977)    in    a   case   study. 

Curing  the  following  24-h  period  ending  0000  GHT  14  Jan- 
uary 1979,  continued  development  is  evident  in  the  low  and 
middle   troposphere.  Central  surface   pressure      (Fig.       10c) 

falls  steadily  tc  977  mb,  which  is  a  decrease  of  18  mb.  If 
one  adopts  the  explosive  criteria  of  Sanders  and  Gyakum 
(1980)  for  a  latitude  of  40  degrees,  the  ston  has  entered 
the  explosive  stage  during  this  period.  The  surface  low 
moves  rapidly  toward  the  northeast,  to  a  position  close  to 
the  maximum  SSI  gradient  (Fig.  11).  Close  proximity  to  this 
gradient  fulfills  another  observation  of  explosive 
development  described  ty  Sanders  and  Gyakum  (1980). 
Analysis  of  the  thickress  contours  (Fig.  10c)  indicates  a 
pattern  of  diffluence  10  degrees  of  longitude  downstream 
from  the  surface  center.  The  diffluence  at  10  degrees 
upstream    is    considerably    greater.  Bjerknes    (1954)       shewed 

that  these  conditions  favor  formation  of  a  diffluent  trough 
in  the  upper  levels  within  the  cyclonic  shear  zone  north  of 
the  jet  streak  axis.  The  NW-SE  tilt  of  the  500  mb  trough 
evident  by  1200  GMT  14  January  (Fig.  16b) ,  combined  with  the 
favorable      superposition      of     the      iet      streak      support      the 
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diffluent  trough  formation  argument.  In  addition,  a  combi- 
nation of  the  superposition  of  the  upper  level  divergence 
zone  over  the  surface  center  and  the  approach  of  the  upper- 
level  thermal  trough  in  clcse  proximity  to  the  500  mb  trough 
promotes  self-amplification  through  the  mid-troposphere. 
Mid-level  vorticity  increases  due  to  the  movement  of  concen- 
trated cyclonic  shear  vorticity  into  the  budget  volume  as 
the  jet  enters  the  diffluent  trough.  As  a  result,  increased 
vertical  development  occurred  from  vorticity  increases  as  a 
result  of    jet    streak    interaction    and   self-amplification. 

Indeed,  the        first     evidence        of        mid-tropospheric 

development  is  visible  in  Fig.  10b.  The  vorticity  pattern 
has  an  expanded  area  of  significant  vorticity  associated 
with   the    mid-tropospheric   circulation.  A   region   of    strong 

vorticity  also  exists  upstream  northeast  of  Japan,  but  that 
center  appears  tc  be  associated  with  the  polar  low  forming 
to  the  west  (Fig.  10c),  as  indicated  earlier.  The  interest- 
ing feature  tc  notice  is  the  absence  of  the  short  wave  that 
was  present  ever  Japan  2U  h  earlier,  which  would  seem  to 
indicate  that  further  explosive  development  from  the  short 
wave   PVA    aloft    is  unlikely. 


39 


At  250  mb  (Fig.  10a)  ,  however,  the  arrival  of  the  jet 
streak  (maximum  velocity  85  m/sec)  in  favorable  superposi- 
tion over  the  surface  system  acts  to  promote  the  lower-level 
circulation,  Diffluence  in  the  region  of  the  jet  stream 
split  and  cyclonic  wind  shear  to  the  left  of  the  jet  axis 
combine  to  induce  positive  circulation  increases  and  strong 
vertical      mcticn     through   the      mid-troposphere.  Increased 

vertical  motion  enhances  vertical  redistributions  of  all 
properties,  including  an  increase  in  the  convergence  in  the 
surface   layer,    and   hence,    explosive   development   continues. 

Inspection  of  satellite  imagery  (Figs.  12  and  13  valid 
1405  GMT  13  January  1979  and  Figs.  14  and  15  valid  at  0023 
GMT      14      January      1979)  shows      a      discernable      comma-shape 

formation.  In  Figs.  12  and  13,  the  southern  boundary  cf  the 
jet  axis  is  delineated  by  subsidence  along  the  cirrus  band 
paralleling  the  lew-level  front.  The  jet  crosses  just  south 
of  the  vortex  at  A.  The  very  ccld  tops  evident  near  the 
vortex  center  and  along  the  frontal  band  indicate  strong 
convection  throughout  the  region.  An  interesting  feature  in 
Figs.  14  and  15  is  the  relatively  low  level  clouds  (labeled 
A)  that  appear  as  fingers  extending  from  the  vortex  center 
and    look      sinilar  to    feeder     bands.         When    coupled      with   the 
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extensive  cirrus  shield  aloft,  the  cloud  pattern  presents  a 
hurricane-like  appearance.  Twc  significant  changes  have 
occurred  during  the  past  10  h.  The  area  of  deep  convection 
has  expanded  markedly  both  near  the  vortex  and  along  the 
frontal  band,  which  indicates  continued  strengthening. 
Second,  the  storm  has  organized  to  the  point  that  clear 
slots  have  formed  in  the  cold  air  behind  the  front,  which 
indicates  the  penetration  of  cold  and  subsiding  air  toward 
the   vortex. 

Curing  the  next  12-h  period  ending  1200  GMT  14  January 
1979,  the  surface  pressure  falls  an  additional  11  mb  to  965 
mb  (Fig.  16c).  Notice  the  shift  of  the  reference  grid  10 
degrees  to  the  ncrth  and  30  degrees  to  the  east  in  this  and 
subsequent      figures.  The      increased      amplitude        of      the 

thickness  ridge  implies  increased  circulation  in  the  lower 
troposphere,  mcst  probably  as  a  result  of  warm  advection 
and  extreme  later.t  heating  from  deep  convection.  Sensible 
heat  input  presumably  weakens  during  the  period  as  the  storm 
center  moves  north  of  the  maximum  SST  gradient  (shewn 
earlier    in   Fig.    1 1)  . 

At  500  mb  (Fig.  16b)  ,  self- amplificat ion  continues  with 
tha      development    of      the    upper-level      trough    and      associated 
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vorticity  maximum.  The  apparent  lack  of  appreciable  vcrtic- 
ity  advection  at  this  level  suggests  that  continued  intensi- 
fication of  the  system  must  be  manifest  in  either  the  lower 
or  upper  levels.  The  mcst  likely  factors  which  serve  to 
maintain  the  deepening  phase  consist  of  a  combination  of 
processes:  jet  stream      interaction   a-      250      mb;      and      the 

extreme  surface  convergence  generating  pockets  of  sustained 
deep   convection. 

An  increase  in  cyclonic  vorticity  appears  at  the  250  mb 
level  during  this  time  frame  (Fig.  16a).  Cyclonic  curvature 
and  shear  to  the  left  of  the  jet  axis  continues  to  promote 
the  circulation  discussed  above.  The  wind  speeds  in  the  jet 
have  diminished  slightly  in  intensity  due  to  the  decreased 
baroclinicity   through  the    lower    troposphere. 

Strong  development  of  the  system  during  the  previous  12 
h  is  indicated  in  Figs.  17  and  18  (valid  at  1205  SHT  14  Jan- 
uary 1979).  Vertical  velocities  within  the  storm  volume 
increase  during  the  period,  as  deep  convection  becomes 
intense  and  widespread  surrounding  the  vortex  and  along  the 
entire  frontal  band.  The  growth  of  convective  activity 
indicates  that  development  continues  upward  through  the 
troposphere      during    the      period,         as      the    jet      crosses      the 


42 


frontal  band  at  a  more  northerly  point.  In  addition,  the 
"hook"  or  tip  of  the  clcud  band  in  the  northwest  sector 
develops  markedly.  In  typical  cases  (Weldon,  1977)  ,  these 
conditions    indicate    the    onset   of    the   mature   stage. 

Interestingly,  the  next  6-h  period  ending  1800  GMT  14 
January  1979  marks  the  most  explosive  period  of  the  storm's 
development  cycle.  The  surface  pressure  (Fig.  19c)  falls  18 
mb       (averaged        between      FGGE      and      NMC        analyses).  Ship 

observations  at  1200  and  1800  GMT  were  limited  in  bcth 
analyses.  It  is  possible  that  seme  of  the  extreme  deepening 
observed  during  this  6-h  period  actually  occurred  before 
1200  GMT,  bur  there  were  inadequate  ship  observations  to 
detect  it.  Thus,  the  magnitude  of  the  surface  pressure  fall 
may  te  less  than  mentioned  above  or  presented  in  Fig.  3b. 
In  any  event,  the  nature  of  the  explosive  deepening  has 
caused  the  thickness  ridge  to  double  in  amplitude  alone  an 
east-west  orientation,  indicating  that  cold  air  has  begun  to 
wrap   around      the     vortex      at      the      surface.  In     the      mid- 

troposphere  (Fig.  19b),  the  positive  vorticity  center  elon- 
gates along  the  flew  and  intensifies  slightly.  The 
elongation  of  the  vorticity  center  and  the  continued  digging 
of    the   mid-lev€l   trough    indicate    the   system   should   weaken    in 
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the  near  future  (Weldon,  1977).  At  250  mb  (Pig.  19a) ,  evi- 
dence still  suggests  that  the  jet  stream  is  in  favorable 
superposition  over  the  surface  system  and  is  continuing  to 
induce  vertical  lotion  and  associated  low-level  convergence. 
Although  the  jet  remains  fairly  stationary  and  strong,  the 
surface  lew  continues  northeast  at  a  high  speed,  which 
suggests    decreased    jet   influence    in    future   development. 

A  look  at  the  satellite  imagery  (Figs.  20  and  21)  valid 
at  2124  GMT  14  January  1979  shows  the  storm  at  maximum 
development.  The  dry  tcngue  aloft  has  penetrated  to  the 
vortex  center.  A  solid  deck  of  st ratocumulus  to  the  south- 
east of  the  vortex  center  exists  in  the  strong  cyclonic 
surface  flow.  An  interesting  "hole"  appears  to  the  east  cf 
the  vertex  in  the  dense  frontal  band  clouds,  which  is  not 
due  to  orographic  effects  as  the  area  is  still  south  of  the 
Aleutian  chain.  A  detailed  look  at  the  preceding  satellite 
picture  (Fig.  17a)  shows  a  clear  slot  formed  in  about  the 
same  relative  position.  Thus,  the  feature  is  associated 
with  the  rotation  of  the  cloud  band  around  this  dry  tongue 
in  the  strcng  cyclonic  flew  at  the  surface,  and  gives  the 
appearance  of  the  formation  of  a  secondary  vortex  in  the 
frontal    band.         Further    inspection   of   this    imagery    shows   the 
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continuance  of  strong  vertical  motion  induced  by  dee?  con- 
vection alcng  the  frontal  tand,  and  to  the  north  and  east  of 
the  vortex  center.  The  penetration  of  cold  air  into  the 
vortex  is  visible  and  signals  the  start  of  the  occlusion 
phase. 

At  0000  GMT  15  January  1979,  the  surface  circulation  is 
dominated  by  the  surface  lew  over  the  Aleutian  Islands  (Fig. 
22c)       with      a   central      pressure    at      947   mb.  The   explosive 

period  has  ended  with  the  onset  of  tha  occlusion  process. 
The  thickness  gradient  has  diminished  markedly,  and  a  polar 
low  is  organized  in  the  cold  air  to  the  west,  while  a  new 
disturbance  forms  alcng  the  baroclinic  zone  to  the  south- 
west. At  500  mb  (Fig.  22b)  ,  a  closed  circulation  has 
formed,  which  indicates  cold  air  influx  through  most  of  the 
troposphere.  A  bread  positive  vcrticity  pattern  exists,  but 
negative  vorticity  advection  into  the  region  is  inducing 
subsidence  and  filling.  The  system  is  vertically  stacked, 
which  also  indicates  the  fully  mature  stage.  At  250  mb 
(Fig.  22a),  development  cf  the  trough  is  evident  from  the 
cyclonic  rotation  of  the  wind  vectors.  The  jet  has  dimin- 
ished in  intensity  due  to  the  reduced  low-level  baroclmic- 
ity  and  the  migration  of  the  surface  system  north  of  the  jet 
axis. 
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During  the  next  24  hcurs  ending  0000  GMT  16  January 
1979,  the  occlusion  process  continues.  The  surface  system 
(Fig.  23c)  has  stalled  over  the  west  Alaskan  coast  and 
filled  to  955  mb .  Ship  reports  in  the  vicinity  show  35  ft 
seas  with  30  ft  swell  and  indicate  that  the  storm  is  still  a 
potent  force.  Other  signs  indicate,  however,  that  the  storm 
is  beginning  to  dissipate.  The  decrease  in  the  thickness 
gradient  indicates  that  the  lower  tropospheric  front  is 
occluded.  The  500  mb  absolute  vorticity  (Fig.  23b)  decrease 
indicates  a  large  reduction  in  relative  vorticity  from  neg- 
ative vorticity  advection  and  the  discontinuance  of  the  jet 
stream  influence  (Fig.  23a).  Extreme  height  falls  are 
reflected  in  the  250  mb  flow  as  the  wind  field  becomes 
closed,  indicating  an  upward  cyclonic  development  through 
this   level. 

Inspection  cf  the  satellite  imagery  for  this  period 
(Figs.  24  and  25)  valid  at  1005  GMT  15  January  1979  shows  a 
fully  mature  occluded  system.  The  penetration  of  cold  air 
to  the  vortex  center  is  evident  along  with  a  marked  decrease 
in  deep  convective  activity  near  the  vortex  and  along  tha 
frontal  band,  except  in  the  northern  fringe  where  the  strong 
cyclonic    flow   is    still      producing  some   low-level   convergence 
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and  associated  deep  convection.  The  developing  low  to  the 
south   is    just    visible   in    the   imagery. 

By  0000  GMT  17  January  1979,  the  surface  system  has 
filled  to  962  mfc  and  the  thickness  gradient  is  small  (Fig. 
26c)  .  At  500  mb  (Fig.  26b)  ,  the  low  has  combined  with  a 
stationary  feature  to  the  west  to  produce  a  vast  area  of  low 
heights.  Positive  vorticity  has  decreased  as  the  vertical 
gradients  have  deteriorated  markedly  while  the  system  slowly 
dies.  At      250     mb      (Fig.        26a),        the      increased      closed 

circulation  indicates  vertical  stacking  through  the  whole 
troposphere.  These   features      suggest  that      the   storm      has 

reached  the  dissipating  phase,  which  marks  the  termination 
of   the   syncptic    pcrtion    of   the  analysis. 

It  is  interesting  to  notice  that  during  the  period  of 
lowest  central  fressure  (0000  GMT  15  January  1979),  ship 
reports  in  the  region  indicate  4  0  ft  waves,  20  ft  swell  and 
38  kt  winds.  The  dramatic  deepening  of  this  and  similar 
storms  is  net  without  its  conse guences.  These  storms  pose 
significant  threats  to  ships  and  land  facilities  withic- 
their  paths.  The  intent  of  the  next  chapter  is  to  pursue 
those  processes  gualit at ively  described  in  this  chapter  by 
mass    and      vorticity    budget  analyses.        The      budget   framework 
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provides      a      computational   base      through      which      significant 
aspects   of    the   storm   are    explored. 
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V.        MASS    BUDGET    ANALYSIS 

A.        GENERAL 

It  is  clear  from  the  synoptic  discussion  (Chapter  4) 
that  the  key  factors  nec€ssary  for  cyclone  development  are 
inherently  related  tc  dynamical  processes  occurring  in  the 
lower  tropospheric  regime.  One  such  vital  process  is  the 
destabilizaticn  cf  the  surface  layer.  Another  includes  the 
horizontal  and  vertical  mass  circulation  structure,  as  pres- 
sure falls  and  resultant  storm  intensity  are  directly  depen- 
dent upon  the  nature  and  magnitude  of  the  mass  flux  into  and 
out  cf  the  budget  volume. 

The  mass  budget  pcrticn  of  this  chapter  includes  a  quan- 
titative approach  to  determine  the  destabilizat ion  and  hori- 
zontal and  vertical  mass  transport  within  the  budget  volume. 
The  stability  analysis  is  based  on  area-averaged  potential 
temperatures  fcr  the  layer.  For  example,  the  layer  labelled 
775  mb  represents  the  850-700  mb  layer.  The  lateral  mass 
transport  is  also  calculated  for  the  layer  and  is  averaged 
over  time  periods.  Vertical  mass  transport  is  inferred  from 
the  vertical  velocity  (omega)  fields.   Time  periods  refer  to 
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the  12-h  tine  increment  between  synoptic  times.  For  exam- 
ple, period  1206  refers  to  0000-1200  GMT.  Radial  distances 
of  6  and  10  degrees  of  latitude  specify  the  horizontal 
budget  area  and  correspond  to  360  and  600  n  mi  radial  dis- 
tances from  storm  center.  These  particular  radii  were  cho- 
sen after  inspection  of  radial  sections  depicting  the 
storm's  major  circulation  features.  The  inner  radius  (6) 
incorporates  the  zone  of  maximum  cyclone  circulation  associ- 
ated with  the  strong  two-layer  inflow/outflow  regime  (subse- 
quently addressed).  The  outer  radius  (10),  in  contrast, 
shows   the    weaker    storm  circulation   near  its    periphery. 

B.        POTENTIAL    TEMPERATURE    AND    STATIC    STABILITY    ANALYSIS 

The  steady  flow  of  polar  continental  air  from  -he  Sibe- 
rian mainland  over  the  warm  waters  of  the  western  Pacific 
was  described  in  Chapter  4.  Des tabilization  of  the  boundary 
layer  in  this  manner  is  a  critical  factor  for  oceanic  cyclo- 
genesis,  as  demonstrated  by  Roman  (198  1)  and  Sandgathe 
(1981).  The  amount  of  dest abi  lizat ion  is  contingent  upon 
the  latent  heat  release  in  the  layer  and  the  sensible  heat 
flux  from  the  ocean.  These  quantities  are  dependent  to  a 
large  degree  on  the  air-sea  temperature  difference,  and  the 
magnitude    of   the     lapse    rate   through   the      lower   troposphere. 
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Roman  (1981)  found  stability  decreases  of  nearly  2  K  /  100 
mb  fully  12  h  before  the  onset  of  maximum  deepening  using  a 
budget      volume   analysis.  Numerical      modal    simulations      by 

Sandgathe  (1981)  showed  a  similar  trend,  with  magnitudes 
roughly  half  that  observed  by  Roman  (198  1)  for  a  western 
Atlantic  storm.  An  analysis  of  the  destabilization  of  the 
lower  troposphere  in  this  case  will  be  described  using  the 
potential    temperature   vertical-time    sections. 

Inspection      of        the      time      sections        of      area-averaged 
potential    temperature   for      radii    6    and   10     (Fig.         27a   and   b) 
shows  similar   interesting    features.  Notice   that   potential 

temperature  values  in  this  figure  are  presented  for  12-h 
synoptic  times.  Lower  trcpospheric  (1000-750  mb)  and  upper 
tropcspheric  (450-175  mb)  potential  temperature  changes 
imply  warming  during  the  early  time  periods.  In  the  mid- 
tropcsphare  (750-450  mb)  ,  potential  temperature  changes  lit- 
tle in  the  early  time  periods.  Of  interest  is  the  apparent 
lack  of  a  well-defined  trcpopause  (indicated  by  a  dense 
vertical  packing  of  isent ropes)  during  tha  early  time  peri- 
ods. Possible  factors  accounting  for  this  observation 
include  the  iirpact  of  jet  streak  propagation  (described  in 
Chapter    4)       into      the    budget    volume,      or      the    characteristic 
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upward  slcpe  cf  the  tropopause  toward  the  tropics,  cr  a  com- 
bination of  both.  At  later  time  periods,  the  entire  budget 
volume  is  seen  to  cool,  with  the  most  significant  decrease 
occurring  below  350  mb,  as  the  cyclone  takes  on  familiar 
extratropical      ccld-core    characteristics.  The      tropopause 

becomes  well  established  by  0000  GMT  15  January,  which 
coincides  with  the  movement  of  the  budget  volume  cut  of  the 
region  of  the  taroclinic  zone  and  into  the  colder  polar  air. 
Several  substantial  potential  temperature  fluctuations  exist 
during  the  1800  GMT  times  (not  shown)  and  are  addressed 
later. 

To  convert  potential  temperature  changes  to  static  sta- 
bility features,  a  procedure  similar  to  that  employed  by 
Sandgathe  (1981)  will  be  used.  He  defines  the  static  sta- 
bility of  the  layer  as  directly  proportional  to  the  specific 
volume  of  the  layer  and  tc  the  change  in  potential  tempera- 
ture between  pressure  surfaces  (d9/dP)  ,  and  inversely  pro- 
portional to  the  average  potential  temperature  of  the  layer. 
For  periods  cf  interest  here,  the  static  stability  of  the 
1000-500  mb  layer  will  be  explored.  The  specific  volume  of 
the  layer  is  assumed  constant  with  time  (variations  are  less 
than    3%) .       The   potential    temperature    fields    (Fig.      27)     imply 
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rather  uniform  average  temperature  structure  for  the  lower 
troposphere  through  time.  Thus,  the  static  stability  varia- 
tions are  predominantly  a  function  of  the  change  in  the 
lapse  rate   with   time. 

The  static  stability  trends  for  inner  and  outer  radii 
are  shown  in  Fig.  28a.  Destabi lization  of  the  layer  occurs 
for  both  radii  until  1200  GMT  14  January.  The  increase  in 
stability  after  1200  GMT  14  January  is  a  reflection  of  cold 
advection  in  the  lower  troposphere.  The  marked  decrease  in 
stability  during  the  early  periods  (Fig.  28a)  is  associated 
with  an  increase  in  the  surface  layer  potential  temperature. 
Thus,  surface  instability  promotes  the  generation  of  organ- 
ized convection  in  the  storm  region,  and  serves  to  foster 
the  cyclonic  circulation  tendency.  k  similar  trend  was 
observed  by  Sandgathe  (1981)  in  his  numerical  simulations 
(Fig.      28b).  Direct   comparison   cf   observed      and   simulated 

stability  magnitudes  is  not  feasible  due  to  the  significant 
differences  in  the  horizontal  scales  over  whica  the  data  was 
averaged.  It  can  be  seen,  however,  that  the  trends  in  the 
observaticnal  data  match  the  simulated  data  very  closely, 
and  show  the  characteristics  of  one  of  the  main  ingredients 
that   is   favorable  for  oceanic   cy clogenesis . 
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One  interesting  feature  of  Fig.  27a  deserves  special 
scrutiny.  The  dramatic  fluctuations  in  potential  tempera- 
ture (8  K)  between  1200  SHI  14  January  and  0000  GMT  15  Janu- 
ary between  800  and  500  mt  is  remarkable  for  such  a  small 
time    interval.  In    general,      the   atmosphere      requires    much 

longer  periods  for  such  largs  fluctuations.  Inspection  of 
85  0-500  mb  potential  temperature  fields  for  both  radii  for 
times  immediately  preceding  and  subsequent  to  1200  GMT  shows 
significant  warming  evident  in  the  1800  GMT  field.  In  addi- 
tion, potential  temperature  biases  of  smaller  magnitudes 
were  evident  at  the  1800  and  0600  GMT  times  on  other  days 
(Fig.  27a  and  b) .  The  nature  of  these  (1800  and  0600  GMT) 
potential  temperature  fluctuations  suggests  a  problem  with 
the  FGGE  data  assimilation  schemes,  and  its  impact  on  this 
study  is    specifically  addressed    in    Appendix   A. 

C.        HORIZONTAL    MASS     FLUX 

The  system  horizontal  mass  transport  analysis  consists 
of  vertical  time  sections  for  radii  6  and  10  (Fig.  29a  and 
b)  in  terms  of  g/sec  for  the  100  mb  layer.  Flux  values  are 
corrected  for  mean  cyclone  movement  and  residual  imbalances 
(see  Chapter  3).  Time  periods  refer  to  the  12-h  period 
between  synoptic  times.  For  example,  period  1206  refers  to 
0000-1200    GMT    12   January. 
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Inspection  of  Fig.  29a  for  radius  six  shows  a  two-layer 
regime  consisting  of  strong  low-level  mass  influx  associated 
with  surface  layer  convergence  below  a  layer  of  mass  outflow 
associated      with    middle      and     upper-level    divergence.  The 

strength  cf  the  mass  flew  is  seen  to  increase  in  both 
regimes  as  the  storm  develops,  with  a  distinct  inflow  maxi- 
mum near  925  mb  and  an  outflow  maximum  near  359  tnb.  The 
inflcw  maximum  (3.27  X  10  *  *  13)  at  1313  GMT  is  roughly 
twice  the  cutflcw  maximum  (1.56  X  10  *  *  13).  Although  the 
lcw^level  inflow  is  intense,  it  is  restricted  to  a  rela- 
tively narrow  layer.  The  deeper,  though  weaker,  upper-level 
outflow  causes  a  net  vertically  integrated  mass  loss  within 
the  column.  The  net  hori2cntal  mass  outflow  in  the  vertical 
columns  (Fig.  29a)  is  reflected  in  the  surface  pressure 
falls    (Fig.    2b)    during   the  developing    periods. 

The  temporal  variation  in  the  level  of  non-divergence 
(LND)  and  the  associated  inflow  maxima  (Fig.  29a)  between 
1306  and  1406  SHT  indicates  vertical  growth  of  the 
convergent   layer     as    intense   cyclogenesis    occurs.  3y    1306 

GMT,  when  closed  isobars  become  evident  at  the  surface,  the 
convergent  layer  begins  tc  grow  as  the  surface  circulation 
is    established.        The   growth   of      the   convergent    layer   during 
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cyclcgenesis  is  consistent  with  results  found  by  Roman 
(1981).  The  radius  6  LND  reaches  maximum  height  (at  1406 
GMT)  soon  after  the  occurrence  of  the  maximum  surface  layer 
transport  and  just  before  the  period  of  lowest  surface  pres- 
sure (1500  GMT),  and  then  exhibits  a  gradual  reduction  with 
time   as   the   storm  begins    tc  dissipate. 

Comparable  features  exist  within  radius  10  (Fig.  29b) 
but  interesting  differences  exist.  The  mass  buildup  in  the 
outer  volume  (1318  GMT)  shows  a  greater  relative  increase, 
which  indicates  differential  transport  between  the  two 
radii,  since  the  inner  volume  increases  are  incorporated  in 
the  cuter  volume  tendency.  The  differential  mass  transport 
can  te  explained  by  the  sarlier  influence  of  the  migratory 
jet  streak  (discussed  in  Chapter  4)  as  it  propagates  first 
into      the      cuter     budget      volume.  The      decrease      in      mass 

transport  after  1318  GMT  reflects  the  decrease  in  jet  streak 
influence  as  it  propagates  through  the  inner  volume  and  dis- 
sipates (see  Chapter  4).  These  patterns  substantiate  the 
intensity  of  assccLatsd  "weather"  observed  by  stations  hun- 
dreds of  n  mi  from  the  stcrm  center  and  clearly  illustrate 
the  immense  scale  of  these  meteorological  events.  It  should 
also    be   noticed      that   the    increase    in      low-level   circulation 
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during  the  last  few  time  periods  at  radius  10  (Fig.  29b)  is 
attributed  tc  the  incorporation  of  the  developing 
circulation   to   the  south     (see   Chapter   4). 

Several  interesting  features  of  the  development  stage  of 
this  cyclone  can  be  observed  by  inspection  of  -he  radial- 
vertical  sections  of  the  corrected  lateral  mass  flux.  A 
sequence  cf  fcur  times  was  chosen  (0000  and  1200  GMT  14  Jan- 
uary and  0000  and  1200  GMT  15  January)  illustrated  in  Figs. 
30    and  31    respectively. 

At  the  outer  radii  in  Fig.  30a,  the  mass  transport  is  a 
maximum  near  radius  11.  The  LND  is  relatively  constant  in 
height  through  all  radii.  A  broad  region  of  strong  upper 
level  divergence  is  also  evident  near  the  outer  radius. 
Notice  the  significant  change  twelve  hours  later  (Fig. 
30b) .  An  additional  layer  of  surface  convergence  near  775 
mb  approaches  the  storm  center  from  the  outer  radii.  This 
feature  suggests  that  a  major  significant  change  in  the 
circulation  has  occurred  during  the  past  twelve  hoars.  Nor- 
mal cyclone  development  might  explain  a  portion  cf  the 
increase,  but  the  explosive  development  indicated  during 
this  time  frcm  the  rapid  pressure  falls  suggests  seme 
additional    forcing   present  in   the   budget   volume. 
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In  the  subsequent  12-h  period  (Fig.  31a) ,  the  surface 
transport  maximum  decreases  and  moves  well  inside  radius  6. 
Divergence  aloft  is  well  established  and  the  positive 
transport  inside  radius  3  between  400-125  mb  suggests  net 
subsidence  and  adiabatic  warming  near  the  cyclone  core, 
which  arrests  at  least  to  some  small  degree  the  cold-cere 
formation  process.  Transport  values  markedly  decrease  12  h 
later    (Fig.    31b)    and    suggest    cyclone   dissipation. 

These  circulation  features  may  be  compared  with  major 
characteristics  cf  tropical  storm  intensification  (e.g.  Hol- 
land and  Merrill,  1983).  Explosive  development  in  both 
cases  appears  to  be  characterized  by  large  circulation  and 
mass  transport  increases  at  the  outer  radii,  with  subsequent 
inner  radial  development  after  organized  transport  has  moved 
near  the  inner  region  with  deep  convection.  The  inner 
radial  development  is  also  largely  dependent  on  the 
formation  of  a  broad  and  strong  outflow  regime  in  the  upper 
troposphere.  As  these  significant  features  form,  organized 
convection  dominates  the  inner  radii  and  triggers  the  rapid 
pressure  falls  due  to  large  latent  heat  inputs  into  the 
lower  troposphere.  Large  pressure  falls  set  the  stage  for 
increased      winds,      circulation      etc..  The   storm      rapidly 
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develops  with  warming  near  the  central  core  until  the 
surface  system  becomes  detached  from  the  upper-level  venti- 
lation and  divergence  source.  The  significant  observed  dif- 
ference in  the  formation  processes  lies  in  the  expanded 
radial  scale  of  the  extratropical  system  compared  to  the 
mesoscale      structure    of      the     tropical   system.  Otherwise, 

upper-level  ventilation  features,  inner  radius  organized 
deep  convection,  weak  stability  of  the  surface  layer  and 
large  outer  region  transports  in  the  early  development  peri- 
ods suggest  reasonable  similarities.  The  conclusions  cited 
above  are  cnly  suggestions  from  the  available  FGGE  fields 
and  will  not  be  explored  in  detail,  as  these  aspects  are  net 
considered   a    primary    objective  of   this   study. 

D.        VERTICAL    MASS   TRANSPORT 

The  storm  vertical  mass  transport  analysis  consists  of 
vertical  time  sections  of  omega  fields  (Figs.  32,  33  and  34) 
for  radius  6  and  10.  Three  distinct  omega  fields  are  calcu- 
lated by  slightly  different  methods  and  the  differences  are 
examined.  The  omega  fields  discussed  below  represent  the 
vertical  mass  transport  and  necessarily  reflect  periods  of 
intense  horizontal  mass  transport  through  continuity  (see 
Chapter    3) . 
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The  differences  between  the  three  omega  fields  are  shown 
clearly   in   Figs.        32,      33     and    34.  The    budget   calculated 

(QLD)  omega  fields  (Fig.  32)  are  area-averaged  values 
derived  from  the  lateral  transport  through  mass  continuity. 
The  FGGE  Level  Ill-b  vertical  motion  estimates  (Fig.  33)  are 
derived  as  pcint  estimates  without  diabatic  internal  adjust- 
ments. Omega  fields  calculated  at  grid  points  using  the 
kinematic  method,  which  implicitly  includes  small-scale 
diabatic    effects,   are  shown  in    (Fig.    34). 

The  difference  in  magnitude  and  level  of  the  vertical 
moticn  maxima  between  the  three  vertical  motion  fields  is 
interesting.  The  QLD  (Fig.  32)  and  FGGE  (Fig.  33)  vertical 
velocity  maxima  match  fairly  closely  in  time  of  occurrence 
and  magnitude,  but  differ  considerably  in  level.  Compensa- 
tion for  the  residual  imbalances  in  the  budget  calculations 
(discussed  in  Chapter  3)  are  applied  layer  by  layer  (highest 
layers  are  weighted  more)  in  the  QLD  omega  field  derivation. 
Thus,  the  QLD  omega  magnitudes  are  comparable  but  the  level 
of   the   maxima    is    biased    toward  the    higher    layers    (Fig.    32). 

The  kinematic  vertical  velocities  differ  both  in  magni- 
tude and  level.  This  is  not  surprising  considering  the  val- 
ues   include      small-scale    diabatic      effects.         Paegle       (1983) 
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demonstrated  similar  results  and  found  that  the  FGGE  derived 
fields  were  rcughly  1/3  the  magnitude  of  the  kinemat ically 
derived  fields,  which  is  consistent  with  a  comparison  of 
Figs.  33  acd  34  (the  ratio  is  closer  to  2/3  here).  For  the 
purpose  of  this  study,  the  kinematic  omega  fields  are 
considered  to  best  represent  the  actual  vertical  motion. 

A  comparison  of  the  vertical  velocity  fields  with  the 
horizontal  mass  transport  shows  consistent  features.  The 
vertical  velocity  maxima  (Fig.  34a)  correlate  well  tempo- 
rally with  the  surface  convergence  and  upper  level 
divergence  maxima  (Fig.  29a)  ,  and  with  periods  of  rapid 
surface  deepening  (Fig.  2b).  In  addition,  the  more  intense 
vertical  motion  at  radius  6  (Figs.  32,  33  and  34),  which  is 
roughly  1.5  times  greater  than  at  radius  10,  agrees  well 
with  the  satellite  discussion  which  depicted  extreme  amounts 
of  deep  convection  in  the  inner  radius  (see  Chapter  4) . 

E.   CONCLUSIONS 

In  summary,  stability  and  QLD  mass  budget  analyses  have 
revealed  the  following  insights: 

•  Surface  stability  decreases  during  the  early  time  oeri- 
ods  favor  organized  convection; 

•  The  stability  analysis  shows  destabilization  trends 
consistent  with  numerical  model  simulations  (Sandgathe, 
1981)  ; 
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•  Evolution  of  cold  core  characteristics  is  evident 
throughout  irost  of  the  troposphere  after  the  surface 
circulation    intensifies; 

•  Maxima  in  surface  layer  convergence  occur  during  peri- 
ods of  intense  vertical  motion  and  strong  upper-level 
divergence,  coincident  with  observed  maximum  cyclone 
deepening   rates     (Fig.    2b); 

•  The  bread  divergent  layer  aloft  contributes  to  central 
pressure  falls  even  though  the  maxima  are  considerably 
smaller  than  the  associated  surface  layer  convergent 
maxima,  indicating  overall  net  mass  outflow  during 
periods   of    storm  intensification; 

•  The  LND  varies  according  to  the  strength  of  the  surface 
convergent  layer,  and  rises  during  the  period  of  explo- 
sive   development; 

•  Vertical  velocity  fields  employing  kinematic  calcula- 
tion methods  represent  well  the  LND  and  vertical  mass 
transport   which      is    a     vital   link     in   the      lifeblood   of 


storm   dynamics. 

Vertical-radial  cross  sections  of  mass  transport  during 
extratrcpica}.  cyclogenesis  suggest  similar  development 
characteristics      to   tropical      storm  rapid      intensifica- 


tion. 

The      next      chapter      seeks      to      explore      the  circulation 
tendencies    as      delineated   by   the    absolute      vorticity   buildup 

within  the    budget  volume    as     the    storm   evolves.  Horizontal 

and   vertical   transport,         and   sources    and    sinks  of  absolute 

vorticity      will        be      explcred      quantitatively  to      further 
delineate    storm   dynamical    features. 
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VI.       VORTICITY    BUDGET    ANALYSIS 

A.       GENERAI 

The  absolute  vorticity  budget  analyses  are  prasented  by 
vertical-time  sections  and  by  selected  vertical  profiles. 
Time  pericds  refer  to  the  12-h  span  between  synoptic  times. 
For  example,  period  1206  refers  to  0000-1200  GMT  12  January. 
In  figures  not  referencing  time  periods,  times  refer  to  the 
specific  time  (GMT)  .  The  775  mb  layer  refers  to  the  layer 
between  the  £50  and  700  mb.  The  sign  convention  is  positive 
for  processes  producing  vcrticity  increases  in  the  budget 
volume  (sources),  and  negative  for  vorticity  decreases 
(sinks)  . 

The  main  feature  of  this  section  includes  a  discussion 
of  the  absolute  vorticity  time  tendency  as  a  measure  of 
storm  development.  Specifically,  an  investigation  is  con- 
ducted intc  the  role  of  horizontal  and  vertical  absolute 
vorticity  transport,  advsction,  generation  and  dissipation 
in  storm  dynamics.  Horizontal  transport  terms  are  separated 
into  mean  and  eddy  mcde  components  to  gain  better  insight 
into   their   specific    contributions      toward    storm    development. 
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Vertical  transport  terms  reflect  vertical  divergence  and 
advection  frcm  cmega  estimates  using  the  kinematic  method. 
Referral  to  Section  B  of  Table  V  is  recommended  as  specific 
budget   terms   are    identified  and    discussed. 

B.        TIME    SECTIONS    ,     TIME    TENDENCIES    AND    SELECTED    PROFILES 

The  dramatic  buildup  cf  absolute  vorticity  in  the  budget 
volume  at  both  radii  over  the  2-day  span  of  intense 
development  is  illustrated  in  Fig.  35.  At  0000  GMT  13  Janu- 
ary (1300  GMT) ,  the  relative  cyclonic  vorticity  in  the  inner 
volume  through  the  600  mb  level  is  not  much  different  from 
the  ccriolis  parameter  (7.82  X  1 0**-5) .  Notice  the  presence 
of  anticyclonic  relative  vorticity  aloft,  which  is  associ- 
ated with  the  weak  and  dissipating  upper  level  ridge.  The 
dramatic  increase  in  absolute  vorticity  by  0000  GMT  15  Janu- 
ary in  the  inner  radius  occurs  from  two  sources.  The  rapid 
movement  to  the  northeast  over  the  period  has  resulted  in  a 
larger  value  cf  planetary  vorticity  (f  =  11.56  X10**-5), 
while  the  increased  storm  circulation  accounts  for  the 
additional   buildup.  Although    the      surface   circulation      is 

strong,  the  naximum  absolute  vorticity  exists  in  the  middle 
troposphere.  Similar  changes  occur  at  radius  10,  although 
the    magnitude    cf    the    absolute    vorticity  increase   is   smaller, 
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and  there  is  a  uniform  increase  through  the  troposphere 
below  the  jet  cere  level.  Specific  factors  contributing  to 
the  vorticity  changes  at  the  different  vertical  levels  will 
be   discussed   as   the    specific   budget   terms    are   analyzed. 

A  different  perspective  of  the  absolute  vorticity 
buildup  in  the  budget  volume  is  presented  in  Fig.  36.  Here, 
vertical  time  sections  of  vorticity  for  both  radii  illus- 
trate several  key  features.  At  radius  6  (Fig.  36a),  a  vor- 
ticity minimuu.  can  be  seen  at  200  mb  during  the  early 
period.  The  increase  with  time  at  lower  levels  is  a  reflec- 
tion of  the  nascent  cyclone  vertical  structure  which  is  com- 
posed of  anticyclonic  vorticity  aloft  over  the  low-level 
developing  system.  After  1212  GMT,  absolute  vorticity 
increases  throughout  the  whole  troposphere  until  1500  GMT. 
Two  significant  features  are  evident  during  these  -cimes.  A 
large  buildup  occurs  in  the  upper  troposphere  (vicinity  225 
mb)  after  1300  GMT,  which  causes  a  -cemporary  reduction  in 
the    vertical    vorticity  gradient.  A   distinct   lower    tropos- 

pheric  maximum  is  evident  at  1500  GMT,  which  occurs  just 
after  the  period  of  intense  circulation  increase,  and  is 
substantiated  by  the  mass  transport  maxima  shown  in  the  mass 
budget      analysis    (Fig.         29a)         and      the   associated      central 
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surface  pressure  falls  illustrated  in  Fig.  2a  of  the  syr.cp- 
tic  discussion.  This  large  maximum  in  vorticity  is  coinci- 
dent with  the  period  of  lowest  central  pressure  and  maximum 
surface  layer  circulation.  Vorticity  values  are  seen  to 
decrease  in  the  middle  and  lower  troposphere  for  24  h  after 
1500  GMT,  but  at  a  significantly  slower  rate  than  the 
buildup.  Values      increase        again        through      the        whole 

troposphere  after  1600  GMT  due  to  the  inclusion  of  the 
developing   system  to    the    scuth. 

Several  important  differences  exist  at  radius  10  (Fig. 
36b).  Vorticity  values  increase  at  a  slower  rate  and  at 
later  times  than  in  the  inner  volume.  This  result,  indicates 
that  although  vorticity  values  increase  in  both  volumes  dur- 
ing the  development  stage,  the  outer  volume  experiences  a 
less  significant  increase  due  to  the  greater  spatial  extent 
and   relatively   smaller   vorticity. 

Absolute  vorticity  time  tendency,  which  is  the  first 
term  in  the  vorticity  budget  eguation  (see  Table  V) ,  is  pre- 
sented in  Fig.  37.  This  display  is  a  time  derivative  of  the 
fialds  depict€d  in  Fig.  36.  Recall  that  time  periods,  when 
used  in  subsequent  figures,  refer  to  the  12-h  time  interval 
between   synoptic   times. 
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A  large  increase  in  the  vorticity  tendency  in  the  inner 
volume  occurs  in  the  upper  levels  from  1218-1318  GMT.  Vor- 
ticity increases  at  a  constant  rate  in  the  surface  layer 
during  the  development  stage.  The  general  decrease  in  ten- 
dency values  through  the  whole  troposphere  after  1506  GMT  is 
consistent  with  the  onset  of  the  occlusion  phase  and  storm 
dissipation.  It  is  interesting  that  a  vorticity  tendency 
minimum  occurs  near  450  tnt  during  a  period  of  large  surface 
pressure  falls  (Fig.  2b)  and  maximum  mass  transport  (Fig. 
29).  No  correlation  to  synoptic  causes  for  this  tendency 
minimum  is  apparent  from  detailed  inspection  of  the  500  mb 
height  and  vorticity  fields  (Figs.  4b,  7b,  10b).  One  inter- 
pretation is  that  a  subsynoptic  scale  system  in  the  wind 
fields  affects  the  vorticity  values  at  these  levels  and 
times.  An  alternative  cause  could  be  linked  to  insufficient 
data   input   into   the    1312    or   1400    GMT   analyses. 

Similar  significant  features  exist  at  radius  10  (Fig. 
37b)  .  The   nature      of    the     buildup      of   absolute      vorticity 

suggests  inherently  similar  mechanisms  are  responsible  for 
the  absolute  vcrticity  buildup  between  inner  and  outer  vol- 
umes. Mechanisms  such  as  the  jet  core  influence,  mid-trc- 
pospheric    self-amplif icaticn,         and   horizontal      and   vertical 
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transport      are      possibilities      explored        in      the      following 
sections. 

C.       LATERAL    TRANSPORT 

The  lateral  transport  of  absolute  vorticity  is  one  of 
the  three  most  significant  forcing  terms  in  the  budget  equa- 
tion (Table  V) .  Transport  into  the  budget  volume  arises 
from  the  mean  convergent  flow  (mean  mode)  and  from 
asymmetries   in   the  cyclonic  flow    (eddy  mode) . 

A  typical  view  of  the  lateral  vorticity  transport  occur- 
ring during  storm  development  is  presented  in  Fig.  38.  The 
curves  are  valid  during  the  maximum  development  period  at 
radius  six.  The  mean  node  and  total  lateral  transport 
curves  are  similar  and  shew  similar  trends.  The  mean  mode 
is  seen  to  dominate  the  total  transport  at  levels  below  800 
mb .  In  the  upper  levels,  the  total  lateral  transport  is 
small,  but  is  the  residual  between  two  strong  oppositely 
contributing  processes.  These  features  constitute  signifi- 
cant differences  from  these  observed  by  Conant  (1982)  in  the 
President's  Day  storm,  where  he  observed  that  the  total  and 
eddy  mode  lateral  transport  curves  were  the  more  strongly 
correlated.  Although  the  curves  presented  in  that  study 
represented   an   earlier   phase   in      the   storm    development,      the 
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general  features  indicate  that  different  mechanisms  serve  -co 
increase  inner  volume  vorticity  levels  within  each  storm 
budget   volume.  A    look    at  the      mean   and    eddy      mode   lateral 

transport  vertical-time  sections  gives  insights  into 
processes    responsible   for    these    differences. 

Another  interesting  feature  is  the  presence  of  a  measur- 
able eddy  mcde  contribution  in  the  lower  troposphere,  which 
is  an  indication  that  processes  other  than  mean  convergence 
are  acting  to  increase  the  vorticity  and  storm  circulation 
within     the   six      degree    vclume.  These      results   should      be 

expected,  since  processes  such  as  intense  convection  in  pre- 
ferred quadrants  of  the  storm,  frontal  asymmetries  and 
elliptical  vortex  structure  may  cause  asymmetric  flow  about 
and  into  the  cyclone.  Thus,  some  eddy  mode  contribution  is 
evident   through   the    whole    troposphere.  Mean    and  eddy    mode 

contributions  to  the  total  vorticity  lateral  transport  are 
discussed    below. 

The  total  lateral  vorticity  transport  time  sections  for 
the  inner  and  outer  budget  volumes  are  presented  in  Fig. 
39.  In  the  inner  volume  (Fig.  39a),  inward  transport  exists 
below  775  mb  through  all  time  periods,  with  alternating 
inward        and      outward        transport      evident        in      the        upper 
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troposphere  (near  350  mb) .  Notice  the  mid-tropospheric  (600 
mb)  outward  transport  (after  15  06  GMT)  as  broad  divergence 
is    established.  The   outward  transport      maximum    at      225    mb 

suggests  strong  divergence  aloft  as  the  jet  stream  dominates 
the  upper  troposphere  flow.  The  surface  layer  regime  is 
particularly  noteworthy  in  that  especially  deep  inward  vor- 
ticity  transport  occurs  during  the  deepening  phase 
(1318-1506  GMT)  through  450  mb  and  contains  a  maximum  (19.2 
X  10**-10/sec**2)  at  925  mb .  This  signature  is  entirely 
different  from  that  observed  by  Conant  (1982)  in  the  Presi- 
dent's Day  storm.  He  showed  a  maximum  at  220  mb  with  rela- 
tively weak  surface  transport.  The  emphasis  on  the  level  of 
maximum  hcrizcntal  lateral  transport  serves  to  focus  on  the 
differences  between  the  underlying  processes  responsible  for 
vorticity  increases  within  the  budgets  of  both  systems.  In 
this  case,  the  lateral  transport  surface  maximum  occurs  in 
temporal  coincidence  with  the  largest  central  pressure  falls 
and  maximum  mass  transport,  as  discussed  in  previous 
chapters. 

The   total      lateral   transport      at   radius      10    (Fig.         39b) 
exhibits      similar     inward      transport   in     the      surface      layer 
although    the    magnitude   is    considerably      less    and   the    maximum 
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occurs  a  time  period  earlier.  Inward  vorticity  transport,  is 
evident  through  the  rest  cf  the  troposphere  (except  abcve 
450  mb  after  1606  GMT)  ,  with  a  secondary  maximum  observed 
between  450  and  175  mb  during  the  1218-1318  GMT  periods. 
The  time  and  level  of  occurrence  of  this  maximum  strongly 
suggests  vcrticity  import  from  jet  streak  interaction.  Sub- 
sequent low-level  inward  transport  occurs  from  circulation 
increases      in   response      tc  this      upper-level    forcing.  The 

earlier  occurrence  of  both  upper  and  lower-level  maxima  in 
the  cuter  volume  suggests  earlier  forcing  from  jet  streak 
interaction  as  it  enters  the  outer  volume  first  and  subse- 
quently migrates  toward  the  inner  volume.  The  increase  in 
the  surface  layer  transport  at  1618  GMT  is  again  due  to  the 
developing   system  to    the    scuth. 

D.       MEAN    MCDE   LATERAL    TRANSPORT 

As  described  in  Chapter  3  (Table  V)  ,  the  absolute  vcr- 
ticity transport  can  be  partitioned  into  mean  and  eddy 
modes.  The  purpose  of  this  partitioning  serves  to  isolate 
contributions  ty  the  mean  flow  (mean  cyclone  convergence  and 
divergence),  and  the  perturbation  flow  (the  covariance  of 
vorticity        and  normal        wind  component        deviations) . 


71 


Specifically,  the  mean  cyclonic  low-level  inflow  and  upper- 
level  outflow  cf  vorticity  is  represented  by  the  mean  mode 
and  reflects  necessarily  the  convergent/divergent  nature  of 
the   lcwer/upper   layers. 

At  radius  six  (Fig.  40a)  ,  the  strong  dual-centered 
transport  maximum,  the  upper-level  outflow  maxima  and  the 
level  of  transition  between  inward  and  outward  transport 
agree  well  bcth  temporally  and  spatially  with  the  mass 
budget  transport  features  discussed  in  Chapter  5.  Periods 
of  maximum  surface  cyclone  deepening  coincide  directly  with 
the  vorticity  transport  maxima  discussed  above.  Differences 
in  mean  mode  transport  between  the  Presidents  Day  storm  and 
the  present  case  arise  due  to  the  significantly  weaker 
surface-layer  lateral  vorticity  transport  in  the  President's 
Day  storm,  whereas  the  upper  layer  transports  are  nearly  the 
same.  Thus,  the  main  difference  in  the  growth  and  vitality 
in  this  storm  is  the  significantly  higher  mean  flow 
transport    of    vorticity   in   the   lower   layers. 

At  radius  10  (Fig.  40b)  ,  absolute  vorticity  lateral 
transport  values  are  relatively  weaker  than  at  radius  6.  A 
single  inward  transport  maximum  in  the  surface  layer  and  a 
coincident    upper-level  outward   maximum      are   evident.         These 
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maxima  occur  earlier,  again,  in  response  to  earlier  dynamic 
forcing  in  the  outer  vclume  by  the  upper-level  eddy 
structure  changes  (subsequently  addressed) .  The  effects  of 
the   new   lew   tc   the  south    are   again    seen  after    1518  GMT. 

E.        EDDY    MCDE    LATERAL   TRANSPORT 

The  eddy  mode  absolute  vorticity  lateral  transport  gives 
insight  into  the  upper  level  vorticity  increases  which  ulti- 
mately affect  the  lew  level  circulation  tendencies.  As 
discussed  above,  the  eddy  mode  component  represents  the  hor- 
izontal transport  due  to  asymmetries  in  the  cyclonic  flow 
induced    by    wind    field  correlated    vorticity   deviations. 

Inspection  of  the  eddy  mode  absolute  vorticity  transport 
for  radius  6  (fig.  4  1a)  shows  considerable  vorticity  import 
in  the  layer  between  750  and  125  mb  after  1206  GMT.  Within 
this  broad  area  of  positive  lateral  vorticity  transport  are 
two  significant  time  periods.  Between  1206-1306  GMT,  a  max- 
imum cf  vorticity  transport  is  seen  to  exist  in  the  450-275 
mb  layer.  Inspection  of  the  3  00  and  400  mb  vorticity  and 
winds  for  1200  and  1212  GMT  shows  weak  ?VA  into  the  budget 
volume,  and  this  is  substantiated  by  the  coincident  maximum 
in  the  calculated  advection  (Fig.  43a)  discussed  subse- 
quently.        These    facts   suggest   vorticity      import   in   the    form 
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of  the  short  wave  that  becomes  established  over  southern 
Japan  during  this  time  frame  (discussed  earlier  in  the 
synoptic  discussion,  Fig.  4b).  The  effect  of  this  perturba- 
tion diminishes  with  time  as  the  cyclone  accelerates  to  the 
northeast  away  from  the  slower  moving  upper-level  trough 
after  1300  GMT,  and  thus,  the  eddy  mode  transport  also 
diminishes.  Vorticity  increases  are  again  evident  after 
1306  GMT  when  the  jet  streak  has  become  well  established  in 
the  budget  volume  (Fig.  7a).  The  influence  of  the  jet 
streak  continues  until  1518  GMT  when  the  jet  core  becomes 
detached  from  the  budget  volume  as  the  surface  center 
migrates    steadily  northward. 

Significant  eddy  mode  contributions  are  restricted  to 
the  upper  levels,  and  are  only  roughly  1/2  those  observed  by 
Conant  (1982)  in  the  President's  Day  storm.  This  is  an 
extremely  valuable  result,  as  Conant  showed  that  the  migra- 
tory short  wave  played  a  major  role  in  that  storm 
development.  In  this  study,  the  migratory  short  wave  is 
shown  to  be  of  minimal  influence,  and  thus  vorticity 
increases  in  the  upper  levels  occur  primarily  from  cyclonic 
shear  vorticity  input  into  the  budget  volume  as  the  jet  cere 
becomes      established.  While      the      eddy      mode      contributes 
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considerably  to  the  initial  storm  development ,  this  effect 
is  net  nearly  as  dominant  during  the  explosive  stages  as  in 
the  President's  Cay  storm.  In  this  case  then,  it  is  logical 
to  conclude  that  explosive  development  occurs  due  to  upper- 
level  jet  forcing  acting  in  concert  with  the  surface 
development  alcng  the  barcclinic  zone.  The  secondary  eddy 
mode  maximum  near  925  mb  at  1418  GMT  is  presumed  to  be  due 
to  asymmetries  in  the  mean  flow  in  specific  quadrants  of  the 
storm  circulation. 

Similar  features  exist  at  radius  10  (Fig.  41b).  As  in 
the  inner  volume,  significant  contributions  occur  at  the  jet 
core  level  with  maxima  cf  similar  magnitude  occurring  at 
coincident  times.  A  weak  maximum  exists  from  1406-1506  GMT 
in  the  surface  layer,  and  as  in  the  6-degrae  volume, 
represents    small   asymmetries   in    the    mean   cyclonic    flow. 

Alternatively,  the  lateral  transport  may  be  expressed  in 
terms  of  lateral  divergence  and  advection  components  through 
the  use  of  vector  identitites  (see  Table  V).  These  terms 
provide  a  more  traditional  view  directly  in  terms  of  the 
vorticity  equation.  The  divergence  component  of  the  parti- 
tioned lateral  vorticity  transport  (Fig.  42a) ,  is  composed 
of    a   two-layer    regime   with   positive    vorticity   production    due 
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to  convergence  in  the  surface  layers  and  negative  vorticity 
production      with    divergence     aloft.  These    features      match 

closely  the  associated  mass  transport  features  discussed  in 
Chapter  5-  The  temporal  and  spatial  occurrence  of  the 
inflow  and  outflow  maxima  agree  well  with  similar  features 
of  the  mean  mode  discussion,  and  with  the  general 
circulation   increases   and    maximum    pressure    falls. 

The  most  interesting  comparison  occurs  between  the  eddy 
mode  and  hcrizcntal  advection  transports  (Figs.  41a  and  b 
and  43a  and  b) .  The  upper  tropospheric  values  match  well 
both  temporally  and  spatially,  as  well  as  in  magnitude. 
However,  the  advection  term  in  the  lower  troposphere  shows  a 
substantial  negative  vorticity  advection  (NVA)  into  the  vol- 
ume below  775  mb,  which  is  quite  different  from  the  negative 
and  positive  eddy  vorticity  transports.  In-depth  examina- 
tion of  the  vorticity  and  wind  fields  in  the  inner  volume 
during  a  period  of  maximum  development  shows  an  asymmetric 
east-west  oriented  vorticity  field  coupled  to  a  strong 
cyclonic  inflow  pattern,  which  causes  a  net  advection  of 
smaller  vorticity  values  into  the  volume  in  the  northeast 
and  southwest  guadrants.  This  net  NVA  is  further  enhanced 
by   the  40-kt   translation    speed  of      the  surface    center    to   the 
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northeast  during  this  period.  The  change  -co  a  westerly, 
more  zonal  wind  component  with  height  accounts  for  the  sign 
reversal  in  the  vorticity  values  above  775  nb.  Notice  that 
the  addition  of  the  advection  and  divergence  terms  is  very 
close  to  the  sum  of  the  mean  and  eddy  mode  lateral 
transports.  These  results  serve  to  verify  the  integrity  of 
the   transport   calculations. 

In  summary,  the  total  horizontal  absolute  vorticity 
transport  has  teen  shown  tc  be  composed  of  a  mean  mode  rep- 
resenting contributions  frcm  the  mean  mass  circulation,  and 
an  eddy  mcde  representing  contributions  from  asymmetries  in 
the  flow  from  advection  and  jet  core  sources.  The  mean  mode 
transport  was  shown  to  be  very  similar  to  the  total 
transport  through  all  levels,  while  the  eddy  mode  transport 
was  shown  to  be  significant  only  above  500  mb.  The  rela- 
tively smaller  magnitude  of  the  eddy  mode  transport  maxima 
compared  to  that  observed  in  the  President's  Day  storm 
suggests  this  term  plays  a  less  significant  rcla  in  the 
storm  evolution. 

F.        VERTICAL    R5DISTRIEUT  ICN 

This  section  investigates  the  role  of  vertical  vorticity 
transport    from      the    upper    and      lower   layers   into      the    middle 
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troposphere  (as  discussed  earlier  in  connection  with  Fig. 
34)  .  In   this      case     the   diffluent      nature      of  the      upward 

evolving  trough  will  be  explained  in  terms  of  this  mid-level 
vorticity  increase.  The  concentration  of  vorticity  in  the 
mid-troposphere  serves  to  increase  mid-level  circulation, 
increase  theriral  advection,  and  amplify  the  trough. 
Vertical  transport  of  vorticity  was  calculated  using  the 
kinematic    omega    fields    (Chapter    5)  , 

The  effect  cf  the  vertical  transport  in  the  vorticity 
budget  is  given  by  the  diveraence  of  the  vertical  transport 
(Fig.  44a  and  b.)  .  This  term  results  from  the  vertical 
derivative  of  the  transport  between  levels  and  represents 
the  redistribution  of  vorticity  upward  through  the 
troposphere.  Notice      the      two-layer     regime      with      strong 

transport  from  the  surface  layer  into  the  middle  and  upper 
troposphere.  The  vertical  transport  of  vorticity  is  a  con- 
sequence of  upper-level  divergence  which  induces  lower-level 
circulation  and  vertical  mction  increases  throughout  the  low 
and  middle  troposphere.  The  maxima  in  both  upper  and  lower 
regimes  are  temporally  coincident  with  periods  of  maximum 
vertical  motion  (see  Chapter  5)  .  Values  in  the  outer  radius 
are  ccnsideratly  smaller  which  is  consistent  with  the  less 
intense    vertical    motions. 
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The      divergent      component       (Fig.  45a)         is      the      main 

contributor  tc  vertical  vcrticity  redistribution.  This  term 
is  simply  the  opposite  of  the  lateral  divergence  term  (Fig. 
42a)  ,  but  differs  slightly  due  to  the  kinematic  vertical 
motion  estimates.  There  is  a  decrease  of  vorticity  below 
the  maximum  in  the  vertical  motion,  and  a  gain  above  this 
maximum.  The  reversal  in  sign  from  the  lateral  divergence 
term  (Fig.  42)  is  consistent  from  continuity,  and  close 
agreement  between  the  temporal  and  spatial  magnitudes  shows 
consistency.  Similar  features  exist  at  radius  10,  although 
the    magnitudes   are  considerably    less   as   expected. 

The  vertical  advection  portion  (Fig.  46a)  plays  a  rela- 
tively minor  role  in  the  tudget .  Positive  vertical  advsc- 
tion  is  seen  tc  occur  between  775  mb  and  350  mb  from 
1318-1506  GMT.  Accumulated  values  of  surface  layer  vortic- 
ity are  advected  by  the  vertical  motion  during  these  times. 
Positive  values  below  and  negative  values  above  450  mb  con- 
tribute to  the  vorticity  buildup  in  the  middle  troposphere 
illustrated  in  Fig.  35.  The  relatively  minor  contribution 
of  vertical  advection  to  the  budget  results  is  consistent 
with    results    found   by   Conant    (1982)  . 
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In  summary,  the  vertical  transport  redistributes  the 
large  vorticity  buildup  in  the  lower  troposphere  into  the 
mid  and  upper  troposphere.  The  divergent  component  was  the 
major  contributor,  as  the  advection  term  was  relatively  weak 
throughout  the  budget  volume.  Processes  at  radius  10  are 
similar    with    considerably    lower    magnitudes. 

G.        SOURCES    AND    SINKS 

Sources  cf  vorticity  cccur  from  the  divergence  and  tilt- 
ing terms,  while  frictional  dissipation  is  a  vorticity  sink. 
Generation  cf  vorticity  by  horizontal  divergence  is  the 
fluid  analog  cf  the  change  in  angular  velocity  of  an  object 
when  it  changes  its  moment  of  inertia  while  conserving  angu- 
lar momentum.  If  horizontal  divergence  occurs  and 
circulation  is  conserved,  the  area  enclosed  by  a  fluid  par- 
cel will  incr€as€  and  the  average  vorticity  of  the  parcel 
must  decrease.  Thus,  the  pattern  of  low-level  convergence 
and  upper-level  divergence  in  a  developing  storm  means  there 
is  a  source  cf  vorticity  at  low  laveis  and  a  sink  in  the 
upper  tropcsphere  (Fig.  4 2a  and  b) .  This  term  is  mathemati- 
cally the  same  as  that  which  appears  in  the  lateral 
transport    partition    (see    Section    B,    Table    V)  . 
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An  important  feature  tc  notice  is  the  relatively  large 
magnitude  of  the  maxima  in  both  regimes  compared  tc  the 
total  lateral  transport.  These  patterns  suggest  -hat  the 
divergence  term  plays  the  dominant  role  in  the  budget.  In 
the  lower  troposphere,  convergence  overwhelms  all  other 
processes,  as  is  expected  from  a  scale  analysis  of  the 
general   vorticity  equation. 

The  tilting  term  (Fig.  47a  and  b)  arises  from  vertical 
vorticity  components  generated  by  the  tilting  of  horizon- 
tally oriented  vorticity  elements  by  a  non-uniform  vertical 
motion  field.  Strong  vertical  shear  and  a  horizontally  var- 
ying omega  field  are  prerequisites  for  contributions  from 
this  source.  It  is  apparent  (from  Fig.  47a  and  b)  that  the 
products  of  vertical  shear  and  the  change  in  the  vertical 
moticn  fields  in  the  horizontal  are  small  compared  to 
lateral   and   vertical    transports. 

Frictional  dissipation  (Fig.  48)  is  assumed  to  occur 
only  in  the  surface  boundary  layer  and  is  parameterized 
using  a  stability  dependent  scheme  (Johnson  and  Downey, 
1976).  The  calculated  values  were  about  1.5  times  greater 
than  those  observed  by  Ccnant  (1982),  as  is  expected  from 
the   larger    intensity    of    the      disturbance    studied   here.         The 
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maximum  at  1506  GMT  has  a  magnitude  very  close  to  the  value 
of  the  low-level  maximum  in  the  total  lateral  transport 
(Fig,  39a).  These  results  suggest  that  frictional  effects 
assume  a  prominact  role  in  the  dissipation  of  vorticity  in 
the  surface  layer.  Values  in  the  outer  volume  (Fig.  48b) 
suggest  minimal  frictional  influences  on  the  vorticity 
averaged   ever   a    larger   region. 

In  summary,  the  divergence  term  was  a  major  factor  in 
the  overall  budget  as  a  scurce  of  cyclonic  vorticity  in  the 
lower  troposphere  and  a  sink  aloft.  The  friction  term  was 
an  important  sink  of  vorticity  in  the  surface  layer,  whereas 
the    tilting   term    had    no    significant   effect   on   the   budget. 

H.        EESIDUALS 

The  residual  in  the  vorticity  budget  contains  the 
effects  of  the  c nit  ted  processes  plus  the  accumulated  errors 
in  the  calculation  of  the  resolved  terms.  Inaccuracies  in 
the  vertical  notion  and  horizontal  wind  fields  contribute  to 
the  physical  residual  components,  whereas  spatial  and  tempo- 
ral finite  differencing  and  grid-point  interpolation  inaccu- 
racies are  causes  for  computational  errors.  A  positive 
vorticity  residual  indicates  an  apparent  vorticity  source  in 
which  the  observed  vorticity  increases  at  a  point  are  larger 
than    is    estimated  by    the    computed   terms. 
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Residuals  for  the  inner  volume  (Fig.  49a)  during  the 
development  period  indicate  generally  small  values  in  the 
middle  and  upper  troposphere,  while  large  vorticity  excesses 
are  evident  at  350  mb  (1306  GMT)  and  at  275  mb  (1406  GMT) . 
Since  these  maxima  occur  during  periods  of  rapid  development 
and  maximum  vertical  motion,  it  seems  likely  that  improper 
representation  cf  the  intense  mesoscale  convection  process 
is  responsible.  That  is,  inaccuracies  in  the  vertical 
moticn  fields  and  therefore  the  vertical  transport  could 
lead  to  underestimation  of  the  vertical  redistribution.  The 
pattern  in  Fig.  49a  is  consistent  with  vorticity  transport 
by  convection,  which  has  been  shewn  to  be  important  in  trop- 
ical circulations  (Shapiro  and  Stevens,  1982).  In  addition, 
residuals  could  also  be  expected  to  be  greater  than  normal 
in  areas  cf  significant  jet  streak  activity  where  the  hori- 
zontal wind  fields  are  generally  the  most  inaccurate.  The 
large  negative  residual  in  the  outer  volume  (Fig.  49b)  is  a 
likely  result  of  this  effect.  Although  computational  resi- 
duals are  difficult  to  isolate  in  an  analysis  such  as  con- 
ducted here,  it  appears  that  contributions  from  this  source 
are  relatively  unimportant  given  the  small  residual  values 
through   most   cf    the    troposphere. 
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I.        CONCLUSIONS 

In  summary,  the  increase  of  absolute  vorticity  in  the 
budget  volume  during  the  development  stage  was  shown  to 
occur      frcm   a      ccmplicated   interaction      of    processes.  The 

period  of  surface  layer  circulation  organization  associated 
with  incipient  development  is  distinctly  different  from  the 
pericd  of  explosive  development  in  which  the  circulation 
increases    throughout    the   troposphere. 

During  the  incipient  stage  (period  1206-1318),  the  lower 
troposphere  circulation  increases  were  shown  to  occur  prima- 
rily from  the  source  divergence  term,  as  the  storm  organized 
along  the    tarcclinic    zone.  Lateral   transport    increased   as 

the  circulation  developed  and  contributed  to  small  net  vor- 
ticity increases  since  the  fricticnal  dissipation  was  also 
enhanced. 

Mid-troposphere  vorticity  values  increased  in  the  earli- 
est period  due  to  weak  lateral  transport  in  the  form  of  PVA 
frcm  the  shcrt  wave  established  ever  Japan  (discussed  in 
detail      in      earlier    sections)  .  This     mid-level      vorticity 

import  acted  to  enhance  the  low-level  organization  in  the 
earliest    stages.  At   later   times      prior    to    1318     GNT,      the 

divergence   of   the  vertical   transport   was    seen    to   account    for 
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a  significant  portion  of  the  mid-level  vorticity  increases. 
An  increase  in  the  deep  divergence  layer  above  the  LND 
served  to    cffset    a   portion  of   the   mid-level   increases. 

In  the  upper  troposphere,  vorticity  values  increased  as 
the  jet  core  propagated  into  the  budget  volume.  Additional 
increases  of  srcaller  magritudes  were  evident  from  vertical 
redistribution  from  the  surface  layer.  These  two  processes 
combined  tc  mere  than  offset  negative  vorticity  due  to  the 
increased      divergence.  Thus,         upper      level      forcing      was 

available    for   surface  development. 

The  most  significant  changes  to  account  for  the  explo- 
sive development  after  1318  GMT  occurred  in  the  upper 
troposphere.  The  small  positive  lateral  transport  was  a 
balance  between  the  mean  mode  (divergence)  and  eddy  racde 
(advection)  contributions.  Large  amounts  of  cyclonic  shear 
vorticity  were  advected  into  the  budget  volume  after  1306 
GMT  as  the  jet  streak  migrated  into  the  budget  volume.  The 
corresponding  divergence  increases  occurred  12  h  later, 
which  served  tc  promote  the  explosive  iow-lavel  development. 
The  maximum  divergence  established  a  time  period  later 
served  to  ventilate  the  system  and  induce  the  maximum  mass 
outflow  aloft,  increasing  the  surface  circulation,  and 
promoting    the   explosive    development. 
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Kid-troposphere  vorticity  increases  after  1318  GMT 
resulted  prinarily  from  the  divergence  of  the  vertical 
transport  and  the  upward  development  of  the  mid-level 
diffluent  trcugh  from  thermal  advection.  The  satellite 
sequence  in  Chapter  4  indicates  that  convection  becomes 
intense  and  widespread  during  the  period  due  to  vertical 
stretching  in  the  mid-trcposphere  induced  by  the  increased 
upper-level   divergence.  Vertical    stretching   acts      no   cool 

the    middle   troposphere  and  destabilize   the    vertical  columns. 

In  the  lower  troposphere,  intense  circulation  increases 
resulted  frcm  the  induced  upper-level  forcing  associated 
with  the  jet  streak  irovenent  into  the  budget  volume.  The 
net  vorticity  increases  results  from  the  excess  of  the 
lateral  transport  over  the  frictional  dissipation.  Although 
the  frictional  effects  increased  at  about  the  same  rate  and 
to  the  same  magnitude  as  the  lateral  transport,  there  was  a 
lag  cf  12  h.  When  the  upper  level  forcing  and  ventilation 
decreased  as  the  surface  center  was  removed  from  the  jet 
influence,    the   system  dissipated    after    1506   GMT. 
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VII.       CONCLUSIONS    AND    RECOMMENDATIONS 

A  synoptic  investigation  was  conducted  of  a  Western 
Pacific  case  of  extratropical  explosive  cyclogenesis  daring 
12-17  January  1979.  Storm  related  mass  and  vorticity  budg- 
ets were  computed  for  the  cyclone  using  the  ECHHF  FGGE 
analyses.  Overall  storm  development  features  were  examined 
and    discussed. 

The   following  significant   results    were   developed: 

•  Maximum      circulation    increases      due      to   rapid      pressure 
falls   match    well  temporally    and      spatially   to   maxima    in 
the    surface     layer    inward   mass      transport    (convergence) 
and    upper-level    outward  transport    (divergence). 

•  PVA  in  the  mid- troposphere  in  the  form  of  short  wave 
interaction,  in  combination  with  thermal  advection 
effects,  enhances  the  organization  of  the  low-level 
circulation    in    the    initial    stages   of    development; 

•  The  initial  short  wave  interaction  decreases  after  the 
initial  period  while  surface  development  continues. 
Explosive      development   occurs      coincidentally   with      the 


Th< 

nui 

de< 

aL_ 

values   increased     after   maximum   storm      development, 

cold    air   advecticn    dominated  the    lower   troposphere. 

Mass  transport  radial  cross  sections  suggest  explosive 
development  characteristics  similar  to  those  observed 
in  tropical  storm  intensification  with  the  followina 
similarities:  1.  lew  static  stability  in  the  surface 
layer,  2.  the  availability  of  an  oceanic  heat  source, 
3.  strong  outflow  (divergence^  aloft,  and  4.  strcng 
convecticn    rear    the    inner  radii; 
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•  The  FGGE  data  set  was  found  to  be  lass  than  optimui  for 
supporting  6-h  QLD  budaet  analyses.  The  0600  and  1800 
GMT  time  periods  contain  spurious  potential  temperature 
biases  affecting  all  of  the  budget  calculations.  The 
problem  also  affects  both  initialized  and  unimtalized 
fields  which  sugaests  inaccurate  analysis  inputs  to  the 
model    first    guess. 

Becommendatiors    for    further    research   include: 

•  The  O'Brien  correction  scheme  applied  to  the  horizontal 
mass  transport  should  be  examined  to  establish  a  mere 
optimum    representation  of  the    vertical   motion    fields; 

•  Further  efforts  should  be  made  to  fully  assess  the 
impact  of  the  FGGE  0600  and  1800  GMT  data  problems  and 
their  effects  on  represented  fields  in  the  FGGE  data 
set. 

•  Energy  diagnostic  budget  calculations  should  be  com- 
pleted on  £his  cyclcne  to  isolate  the  contribution  of 
thermal  advection,  sensible  heat  and  latent  heat  in  the 
explosive  development.  Correct  parameterization  of 
these  inputs  could  serve  as  a  basis  for  further  compar- 
isons between  explosive  extratropical  and  tropical 
storm   formation. 

Accurate    and   reliable  forecasts    of    maritime    explosive    cyclo- 

genesis    will   remain    a      difficult    problem   until    observational 

networks    are    improved,      latent  and    sensible   heat   flux    inputs 

to    storm    dynamics  are   better   parameterized,    and   vertical   and 

horizontal    model   resolutions  are    improved. 
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APPENDIX    A 


FGGE    DATA    EIAS    CHARACTERISTICS 


The  following  information  is  presented  to  document  the 
6-h  FGGE  data  problem  encountered  in  this  study.  As  indi- 
cated in  Chapter  5,  a  significant  high  potential  temperature 
bias  was  observed  in  the  850-500  mb  layer  at  1800  GMT  14 
January.  Several  additional  fluctuations  were  observed  to 
occur  at  0600  and  1800  GMT  on  other  days.  The  nature  of 
these  temperature  biases  radically  affected  the  budget  cal- 
culations and  precluded  the  planned  use  of  a  6-h  interval 
format.  The  nature  of  the  effects  is  shown  in  Figs.  50  and 
51.  The  top  portion  of  Fig.  50  is  a  6-h  representation  of 
the  vorticity  time  tendency  term  at  radius  six,  whereas  the 
bottom  field  is  the  inner  volume  budget  residual.  The  sig- 
nificant 6-h  oscillations  in  both  fields  are  indicative  of 
the  difficulty  encountered  using  QLD  budget  techniques  on 
this    particular   FGGE    data    set. 

An  even  more  striking  example  of  the  complexity  of  the 
6-h  data  problem  is  shown  in  the  inner  volume  stability 
trace  (Fig.  51).  Notice  the  oscillation  between  maxima  and 
minima      indicating        the      alternate        warming      and        cooling 
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tendencies  in  the  6-h  fields.  It  is  obvious  that  a  synoptic 
xime  scale  interpretation  cf  this  figure  is  virtually  impos- 
sible. The  1300-1500  GMT  period  is  especially  interesting. 
The  0600  and  1800  GMT  stability  values  are  markedly  larger 
than  at  0000  and  1200  GM1.  The  increased  stabiitiy  during 
these  times  occurs  from  warm  biases  in  the  mid-level 
potential    temperatures,    as   shown    in   Fig.    27a    and   b. 

To  determine  possible  causes  for  this  6-h  data  problem, 
specific  mid-tropospheric  observational  data  available  as 
input  to  the  FGGE  analysis  scheme  at  1418  GMT  was  studied. 
The  observations  consisted  solely  of  satellite  input  (micro- 
wave and  infrared  wavelengths) .  Close  scrutiny  of  the  sat- 
ellite input  in  comparison  with  the  FGGE  gridpoint  estimates 
revealed        interesting      discrepancies.  The      final        FGGE 

potential  temperature  analysis  was  5-7  K  higher  than  the 
associated  satellite  data  in  the  southwest  quadrant  cf  the 
budget  volume.  In  all  other  regions  the  soundings  and  the 
FGGE  fields  matched  very  closely.  In  addition,  potential 
temperature  biases  of  smaller  magnitudes  were  evident  at  the 
1800    GMT   time   period    on    ether   days    (Fig.    27)  . 

Further  insight  is  gained  by  an  inspection  of  the  unini- 
tialized   height      fields     (through      the   thickness      analyses    in 
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Figs.  16c,  19c  and  22c) .  These  fields  also  show  fluctua- 
tions of  the  same  scale  occurring  at  the  1800  GMT  time 
frame,  as  this  thickness  field  departs  markedly  from  the 
preceding  aud  subsequent  analyses.  It  is  obvious  that  the 
6-h  data  problem  affects  both  initialized  (potential  temper- 
ature) and  uninitialized  (height)  fields.  The  evidsnce  pre- 
sented above  suggests  that  a  problem  existed  in  the  analyzed 
potential  temperature  fields  used  in  the  first  guess,  and 
that  the  nature  cf  the  prcblem  was  not  identified  during  the 
FGGE   data    post-processing    evaluations. 

Further  support  of  these  conclusions  is  presented  by 
Paegle  (1983),  who  summarized  features  of  the  FGGE  data  set 
in  tropical  regions.  Lower  tropospheric  temperature  biases 
on  the  order  cf  2-3  degrees  were  found  for  some  time  peri- 
ods, which  were  presumably  from  the  lack  of  conventional 
data  input  at  the  1800  GMT  times.  In  addition,  the  ECMWF 
forecast  model  was  shewn  tc  be  insensitive  to  diurnal  temp- 
erature changes  over  land.  To  what  degree  this  is  evident 
over   oceanic   areas  is   still  uncertain. 

In  summary,  it  is  believed  that  the  spurious  warming  at 
1800  GMT  between  850  and  500  mb  is  due  to  a  problem  with  the 
FGGE     analysis      scheme        rather      than      a      lack        of      quality 
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observations.  The  inpact  of  this  assessment  on  the  results 
of  this  thesis  cannot  be  fully  ascertained  from  the  prelimi- 
nary cbservaticns  discussed  here.  Detailed  investigation 
into  specific  analysis  fields  used  in  the  first  guess  must 
be  dene  level  by  level  tc  fully  assess  the  impact  of  this 
6-h    data    problem. 
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AFPENDIX    B 
TABLES 


TABLE    I 
Availability   of    FGGE   Data 

(From    Eengsston,    1982) 


PH 
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TABLE    II 
Lis-ir.g   of    Analyses   and    Prognoses    Utilized 

(FGGE=1  ,NMC=2,  FNOC=3) 


LEVEL  TYPE 

sea-level  pressure 


ANAL/PROG 
anal 


ORIGIN 

1,2,3 


sea-level  pressure 


36-h   prog 


2,3 


sea-surface     temperature  anal 


500-mb  heights 


anal 


1,2,3 


9U 


TABLE    III 
Generalized   QLD  Budget    Equation   in   Isobaric   Coordinates 

(After   Wash,     1978) 


100  nb       B       2ir 
F  =        /  /      /        -  fr2  sing  dadS    (-dp) 

1000  nb     0        0 

where  F  is   the  volume   integral   of    the   desired  budget   property   f, 
The  budget   eauation  is 


|f  =  LT(f)   +  VT(f)   +  S(f) 

ot 


where  the  lateral  transport  is 

100  nb  2i 

LT(f)      /  /    -  (U-W)6  frsinS  dcx(-dp)|g 

1000  nb  0 

and  the  vertical  redistribution  is 
100  nb   3   2tt 


/ 


VT(f)  =  '    -  (oof)  r2  sinB  dad6  (-dp) 


1000  nb  0   0 

The  source/sink  tern  is 

100  nb   3   2tt 


/ 


/      I       -  ~   r2  sin3  dadS  (-dp)  . 


1000  nb  0   0 
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TABLE    IV 
QLD    Mass   Budget    Equations  in   Isobaric   Coordinates 

(After  Wash,    1978) 
The  Definition 


M  =       I     -  r2  sinB  dadBdp 

V 
P 


where  f  =  1    . 


The  Budget   Equation 


^  =  LT  +  VT    , 
at 


r.th  r\t~  r\ 


Lateral  Transport 


'T 


=  -   I       ~  (U-W)^  r  sin8s  dadp  |g 


Vertical  Transport 
3R   2tt 

r  /"    i        2 

VT  =    /   /    -  (0)-oO  r  sing  dadS 

J      J  8      B  p 

0   0 
where 


dp.  dpB 

M  =  dT  •  "b  =  IT  • 
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TABLE  V 
QLE  Circulation  Budget  Equations  in  Isobaric  Coordinates 

(After  Wash,  1978) 

Section  A 

100  mb        8        2tt 


/  I    i 


Ca  =  /  /       /         -  ^a  r2  sinB  dctdB    (-dp) 


1000  mb     0       0 


where  C  is  the  absolute  circulation  and  L,      is  absolute  vorticity, 
a  a 


The  budget  equation  is 


-*-a  =  LT(r  )  +  DVT(C  )  +  S(t  ) 
ot        a        a       a 


where   the  lateral   transport  is 
100  mb        3        2tt 

LT(Ca}  /  /       /  g    Q'Vtl   Car   Sin6    dct(~dP)l3 

1000  mb     0        0 

and   the  divergence   of   the  vertical   transport   is 
3        2tt 


•//i 


DVT  (5   )  -     /         I         —  |-  (oo  5   )   r2  sinB   dadB    . 
a  /         J  g  3p  a 


0        0 

The  source/sink   term  is 
100  mb       3        2tt 


s(ca)  =       /  /     /       7ir  r2  sin3  dad^"dP) 


1000  mb     0        0 
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TABLE  V  (Cont.) 


Section  B 


The  partitioned  form  of  the  vorticity  budget  equation  is 


6(C) 
a 

at 


+  LT  (£  )      +      DVT  (5  ) 

a  ,  a 


+    S  (C  ) 

a 


mean  mode 


'       vertical    vertical 
eddy  mode   divergence  advection 


horizontal  horizontal 
divergence  advection 


, 

divergence   tilting  frictional 
term      term    dissipation 


The  above  partitions  make  use  of  Stokes'  theorem 


/£  U  •  m 
a~   . 


■  -   I  |V«?  UdA, 


and  the  divisxon  of  total  flux  (U  £  )  into  divergent  and  advective 

~   a 
components, 


7-C  U  =  5   (7-U)  +  U«V  ?   . 
a  ~    a    ~    ~    a 
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TABLE    VI 
Mear.   and   Eddy    Mode   Transporx  Equations 

(After    Conant,    1982) 


100  mb  2tt 

f  f  1~ ** 

EM(C   )   =  /  /  -  C      (U-W)Q        r  sing  da(-dp) 

a  J  /  g     a     ~   ~  p 

1000  mb  0 


100  mb  2tt 


3 


/   /  i 


MM(£   )   =         /  /  -r    (U-W)^       r  sing  da(-dp)l(. 

a  /  /  g     a     ~    ~   p 

1000  mb        0 


100  mb  2tt 

f  f  1~* *° 

E>I(ga)  /  /  ~~  Sa    ?"?)3        r  sinS   dcx(-dP) 

1000  mb  0 


3 


100  mb  2tt 

MM(ga)  /  /  g  S°    (U-W)g        r  sinS   da(-dp)|g 

1000  mb        0 


EM  is  the  eddy  mode  and  MM  is  the  mean  mode. 

2tt 
Note  that   (  )3  =  -jj-   /    (  )  da 

0 


*         — .  a 
and        ()=()-()    is  the  deviation  of  the  property 

from  its  mean  around  a  lateral  boundary. 
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SEPENDIX    C 
FIGURES 
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Figure    1. 


The   Overall   Data   Management    During    FGGE    (From 
Bengsston,    1982). 
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Figure   2. 


Storm   Track  Positions    and   SLP    Analyses.       (A) 
Storm    Positions   in    12-h    Incraaents.       (B)    Sea 
L«vel    Pressure    Analyses    in    12-h    Increments.    Times 
12/0000   and    1200  refer   to   0000    GMT    12   January 
1979. 
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Figure  3. 


36-h    Forecast    SLP   and   Storm   Position   Errors.       (A) 
S-a   Level    Pressure   Forecasts  in    12-h    Increments. 

E)    Position    Errors   Based   on   FGGE    Positions. 

ime    1200    refers  to  0  000   GMT   12    January    1979. 
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Figure  4.      Synoptic    Fields  for  0000   GMT   12    January    1979. 

(A)    winds/Isot  achs    (i/sec).     (B)     500    mb    Absolute 
vcrticity    (DASHED)    X    10**5/sec   and   Heights 
(SCLID)    in   gpm.    (C)    Sea    Level    Pressure    (SOLID)    in 
mb    and   1000-500   mb  Thickness    (DASHED)    in   gpm. 
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Figure  5.  EWSE  Visual  Satellite  Imaaery  for  0041  GMT  12 
January  1979.  Labels  are  jA)  Incipient  System 
(E)  Secondary  DevelcDing  System  (C)  PVA  Area 
upstr< 


earn. 
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Figure  6.      DHSE    Infrared    Satellite    Imagery    for 


January  1979.  labels  are  (A)  Incipient  Svs 
(E)  Secondary  Developing  system  (C)  PVA  Ar 
upstream. 
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Fiaure  7.   Same  as  Figure  4  except  for  0000  GMT  13  January 
1979. 
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Figure   8.       Same    as   Figure    5   except    for   234  1    GMT 
1979.    Label   A    is  Incipient    System. 


12   January 
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Figure  9.       Sams   as   Figure   6  except    for   234  1    GMT    12   January 
1979.    Label  A    is  Incipient    System. 
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Fiqura    10.      Same    as    Figure   4    except    for   0000    GMT    14   January 
1979. 
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Figure   1 1 . 


FNCC    SST    Analysis   Valid    1200. GMT    14   January 
1S79.      "X"    marks   cyclone    position   at   0000    GMT    14 
January.      Isotherms    in    degrees   Celsius. 
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Figure    12. 


Same   as   Figure 
1S79.      Label    A 

Intersection. 


5   except    for    1405    GMT    13   January 
Indicates   Jet-Frontal 
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Figure   13. 


Same    as    Figure   6    except    for    1405    GMT    13   January 
1979.      Label    A   Indicates  Jet-Frontal 
Intersection. 
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Figure    14. 


Same   as    Figure   5   except   for   0023    GMT    14  Januarv 
1S79.      Label    A   Indicates   Feeder    Bands. 
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Figure    15. 


Same   as    Figure  6   except   for  00  23    GMT    14   January 
1979.      Label    A   Indicates    Feeder    Bands. 
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Fiqure    16.      Same    as   Figure  4    except   for    1200    GMT    14   January 
1S79. 
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Figure  17.   Same  as  Figure  5  except  for  1205  GET  14  January 
1979.  J 
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Figure  18. 


Same  as 
1979. 


Figure  6  except  for  1205  GMT  14  January 
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Fiaure  19.   Same  as  Figure  4  except  for  1800  GMT  14  January 
a  1S79. 
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Figure  20.   Same  as  Figure  5  except  for  2124  GMT  14  Jaruary 
1  S79 . 
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Figure  21.   Same  as  Figure  6  except  for  2124  GMT  14  January 
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Fiqure   22.      Same   as    Figure   4   except   for   0000    GMT    15   January 
*  1S79. 
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F^  23.      Same   as   Figure  4    except   for   0000    GMT    16   January 

y  1979. 
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Figure  24. 


Same  as  Figure  5  except  for  1005  GMT  16  Januarv 
1^79. 
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Figure  25.   Same  as  Figure  6  except  for  1005  GMT  16  January 
1S79. 
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Figure 


26         came   as   Figure   4   except   for   0000   GMT    17   January 
1979. 
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Figure    27. 


Area-Averaged    Potential   Temperature   Time 
Sactions     (A)     Radius    6    ( B)     Radius    10.       Values 
in    Degrees   Kelvin.       Time    1200    Refers   to   0000 
12   January. 
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Figure   28 


Stability  Index  (A)  Observed 
Sandgathe,  1981).  Values  in 
100  mb.  Time  1200  Refers  to 
January. 


(B)  Simulated  (From 
Degrees  Kelvin  / 
0  00  0  GMT  12 
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Figure   29. 
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Corrected   Lateral   Mass   Transport    Time    Sections 
(A)    Radius  6     <E)    Radius    10.      Contour    Interval    is 
U0   X    (10    *   *    11).      Dashed    (Negative)    Values 
Indicate    Flux   Cut   of    the   Volume.      Units   are    in 
am/sec-100   mb.      Period    1206    Refers   to    0000-1200 
GMT    12  January. 
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Figure   30-      Radial  Section  of   Lateral    Mass   Transport    (A) 
Time    14  0  0 
in   Figure 


Time    1400    (B)     Time    1412.       Units   and    Contours    as 
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iqure  31.  Radial  Section  of  Lateral  Mass  Transport  (A) 
Time  1500  (B)  Time  1512.  Units  and  Contours 
in  Figure  29. 


s  as 


130 


(fl) 


100 


150- 


200 


cj 


300- 


<o 

CO  400 
500 
700 
850 


SFC 


tf 


'III  %    >    n  N                .*'    N 

i  ;     i  »   \   >■  \        *>    n 

I   '   '  \     *     \  '                         N 
'   '        '  X     \  N 


t       t        f 

'  t 

'  \ 

•  » 

»      » 

/       I 
/        I 

/    \  y 


/      ; 
/ 

/ 
/ 

l 

\ 
\ 
I 

\ 
\ 
\ 
\ 
J 


.;:f:::r:.— .—  — — • 


161B 


1318     1406     M18     1506 

TIME  PERIOD 


1618 


Figure  32.   QIC  Derived  Vertical  Velocities  (A)  Radius  6  (3) 
Radius  10.  Units  are  in  mb/1 000-sec.   Time 
Periods  as  in  figure  29. 
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Figure    33.      FGGE   Derived    Vertical   Velocities    (A)    Radius    6 

(E)    Fadius    10.      Units    are    in    mb/1000-s9c.       Time 
Periods    as   in    Figure    29. 
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Figure   34.      Kir.ematic  Vertical  Velocities     fA)    Radius    6     (B) 
Radius    10.       Urits   are   in   mb/1 000-sec.      Time 
Fericds    as   in    Figure    29. 
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Figure    35.      Budget   Volume    Absolute    Vor-ciciry    for    Radius    6 

and    Badius    10    for    Two   Times.    Time    1300    Z    Refers 
tc    0  000    GMT    13   January. 


134 


SPC 


12m     1212     1300     1312     MOO     1412     1500     1512     1600     1612     !76o 


1700 


Figure  36.      Area-^averaged    Absolute    Vcrticity    Vertical   Time 


Sections    for    (A)    Radius   6    and     (B)    Radius    10. 
Time    1200    Refers   to    0000   GMT    12   January 


are    10**-5/sec. 
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Absolute    Vorticity   Vertical  Time-Tendency 
Sections    for     (A)     Radius    6    and    (B)     Radius    10. 
Uric!    1206    Refers  to   0000-1200    GMT    12    January, 
Units  are   10    *  *   -10/(sec   *.*    2).      Solid 
Ccntcurs    reflect    vorticity   increases. 
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Figure    38. 


Tctal,   Mean    and  Eddy    Mode    Lateral    Vo: 
Transport   During   Maximum    Development 
in   Figure   35.      units   as    in   Figure    36. 
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Figure  39.   Same  as  Figure  37  except  for  Lateral  Vorticiry 
Transport. 
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Figure   40.      Same   as    Figure   37   except    for    Mean   Mode   Lateral 
Transport. 
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Figure  41.   Same  as  Figure  37  except  for  Eddy  Mode  Lateral 
Transport. 
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Figure  42.      Same   as    Figure   37    except    for   Lateral  Vorticity 
Divergence. 
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Figure   43.      Same   as    Figurs   37    axcept    for   Lateral   Vorticity 
Advection. 
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Figure  44. 


Divergence  of   the   Vertical  Transport.      Units, 
times  and  contours  as   m   Figure    37. 
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Figure  45.   Same  as  Figure  37  except  for  Vertical  Vorticity 
Divergence. 


144 


1206 


13i)6  1318 


1106  uia 

TIME  PERIOD 


1518 


1606  1618 


Figure  46.   Same  as  Figure  37  except  for  Vertical  Vorticity 
Advection . 
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Figure  47.   Same  as  Figure  37  except  for  the  Tilting  Term 
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Figure   48.      Same   as    Figure   37   except   for    Frictional 
Dissipation. 
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Figure  49.   Same  as  Figure  37  except  for  Budget.  Residuals. 
Positive  Values  Indicate  Vorticity  Excess. 
Negative  Values  Indicate  Vorticity  Deficits. 
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Figure   50. 


FGGS    6-h    Radius   6    Vorticity   Budget    Fields    for 
(A)    Time    Tendencies     (B)     Buaget    Residuals. 
Period    1203    Refers   tc    0000-0600    GMT    12    January 
and    is  Valid    for   0300   GMT.      Contours   as   in 
Figure   37. 
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Figure   50. 


FGGE    6-h    Radius   6    Vorticity   Budget    Fields    for 
(A)    Time    Tendencies     (B)    Buaget    Residuals. 
Period    1203    Refers   tc   0000-0600    GMT    12    January 
^nd    is   Valid    for    0300    GMT.      Contours   as   in 
Fiaure   37. 
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Figure   51.      FGGE    6-h    Stability   Trace.      Units   and  Times    as   in 
Figure   28. 
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